Quantum simulation of open-system dynamical maps with trapped ions 
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Dynamical maps describe general transformations of the state of a physical system, and their 
iteration can be interpreted as generating a discrete time evolution. Prime examples include clas- 
sical nonlinear systems undergoing transitions to chaos. Quantum mechanical counterparts show 
intriguing phenomena such as dynamical localization on the single particle level. Here we extend the 
concept of dynamical maps to an open-system, many-particle context: We experimentally explore 
the stroboscopic dynamics of a complex many-body spin model by means of a universal quantum 
simulator using up to five ions. In particular, we generate long-range phase coherence of spin by an 
iteration of purely dissipative quantum maps. We also demonstrate the characteristics of competi- 
tion between combined coherent and dissipative non-equilibrium evolution. This opens the door for 
studying many-particle non-equilibrium physics and associated dynamical phase transitions with no 
immediate counterpart in equilibrium condensed matter systems. An error detection and reduction 
toolbox that facilitates the faithful quantum simulation of larger systems is developed as a first step 
in this direction. 

PACS numbers: 32.80.Qk, 64.70.Tg, 37.10. Ty, 03.67.Ac, 42.50.Dv 



I. INTRODUCTION 

Obtaining full control of the dynamics of many-particle 
quantum systems represents a fundamental scientific 
and technological challenge. Impressive experimental 
progress on various physical platforms has been made [1]- 
[12], complemented with the development of a detailed 
quantum control theory |13H16j . Controlling the coherent 
dynamics of systems well-isolated from the environment 
enables, for example, quantum computation in the circuit 
model [17] . But this also allows for digital coherent quan- 
tum simulation with time evolution realized by sequences 
of small Trotter steps [TH] , as demonstrated in recent ex- 
periments pjJl [20] . On the other hand, engineering the 
coupling of a system to its environment - and thus its re- 
sulting dissipative dynamics - introduces new scenarios 
of dissipative quantum state preparation 21 25 , dissipa- 
tive variants of quantum computing and memories (26j- 
128] and non-equilibrium many-body physics [29H31] . Fi- 
nally, complete control of both coherent and dissipative 
dynamics would enable the operation of an open-system 
quantum simulator, whose elementary building blocks 
have been demonstrated recently [25] , and which holds 
the promise to experimentally explore the dynamics of 
novel classes of non-equilibrium multi-particle quantum 
systems. 

The dynamics of these systems is often considered as 
continuous in time, described by many-body Lindblad 
master equations, cf. e.g. [32j . This may be conceived as 
a special instance of a more general setting, where a dis- 



crete time evolution of a system's reduced density matrix 
is generated by concatenated dynamical maps. So far, 
the concept of dynamical maps has proven useful for the 
description of periodically driven classical nonlinear sys- 
tems [55] . and their quantum mechanical counterparts, 
such as the kicked rotor, providing one of the paradig- 
matic models of quantum chaos [34M36j . Remarkable ex- 
periments have been performed with periodically driven 
systems of cold atoms, which have demonstrated some of 
the basic phenomena of quantum chaos such as dynami- 
cal localization [571 - |4"0"] . At present all these studies are 
on the level of single particle physics. 

In Sections II & III below, we present a first experi- 
mental study of open-system many-particle quantum dy- 
namical maps for a complex spin model implemented in 
a linear ion trap quantum computing architecture us- 
ing up to five ions. Building on recent advances in our 
experimental techniques, we harness the building blocks 
of our digital open-system quantum simulator device to 
concatenate elementary quantum maps with high fidelity. 
The controlled open-system maps are enabled by a dig- 
ital simulation strategy, contrasting analog Hamiltonian 
quantum simulation with trapped ions |41H45j . In this 
framework, we observe the key physical phenomena of 
multi-particle dissipative and coherent interactions in the 
stroboscopic system dynamics. In particular, we demon- 
strate the purely dissipative creation of quantum me- 
chanical long-range phase order. Furthermore, we im- 
plement a competition between coherent and dissipative 
many-particle dynamics by alternating sequences of uni- 
tary and non-unitary maps outlined in Fig. [l^,, and ob- 
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serve the destruction of phase coherence as a result. This 
reflects the hallmark feature of a strong coupling non- 
equilibrium phase transition predicted in a closely related 
driven-dissipative model of bosons |46j . 

The presented material is complemented by an exten- 
sive and self-contained collection of appendices, which 
provides further mathematical concepts and technical de- 
tails. 

II. COMPETING DISSIPATIVE AND UNITARY 
DYNAMICS IN A COMPLEX SPIN MODEL 

Two competing, non-commuting contributions to a 
Hamiltonian give rise to a quantum phase transition, if 
the respective ground states of each contribution sepa- 
rately favor states with different symmetries [J? . The 
transition takes place at a critical value of the dimen- 
sionless ratio g of the two competing energy scales. A 
non-equilibrium analog can be achieved in open many- 
body quantum systems, where coherent Hamiltonian and 
dissipative dynamics compete with each other: The role 
of the ground state is played by the stationary state of 
the combined evolution, and the dimensionless ratio g 
is provided by a Hamiltonian energy scale vs. a dissi- 
pative rate. Such a situation has been addressed pre- 
viously theoretically in the context of driven-dissipative 
dynamics of atomic bosons on a lattice [46]: A dissipa- 
tive dynamics can be devised to drive the system from an 
arbitrary initial state with density matrix p- m into a Bosc- 
Einstein condensate with long-range phase coherence as 
the unique, pure "dark" state d of the dissipative evo- 
lution, i.e. pi n — > po = \^P)d (V 1 ! f° r long enough waiting 
time. Supplementing this dynamics with a Hamiltonian 
representing local interactions, being incompatible with 
the dissipative tendency to delocalize the bosons, gives 
rise to a strong coupling dynamical phase transition. It 
shares features of a quantum phase transition in that it 
is driven by the competition of non-commuting quantum 
mechanical operators, and a classical one in that the or- 
dered phase terminates in a strongly mixed state. 

In our experiment, we consider analogous open-system 
dynamics of a quantum spin- 1/2 - or hardcore boson - 
model, realized with trapped ions. A schematic overview 
of the relation of the ionic spin- and the atomic boson 
model is given in Fig. [TJ), whereas a more detailed de- 
scription is provided in Appendices |B|and[D"l The discrete 
time evolution is generated by sequences of dynamical 
or Kraus maps £ w acting on the system's reduced den- 
sity matrix p as illustrated in Fig. [l£i, with time steps 
ti — > tg + i represented by 

K 

p(t e ) ^ p(t t+1 ) = 8®\p(t t )] = J2E^p(h)Ei l) l 

k=l 

The set of Kraus operators {Efp} satisfies 

Ef=i E k )J * E k ] = 1 PU- While the familiar se- 
quences of unitary maps are obtained for a single Kraus 
operator K = 1, dissipative dynamics corresponds to 



multiple Kraus operators K > 1. In particular, the 
continuous time evolution of a Lindblad master equation 
is recovered in the limit of infinitesimal time steps, cf. 
Appendix [SJ The dissipative dynamics studied in our 
spin model is governed by dynamical maps according to 
two-body Kraus operators acting on pairs of neighboring 
spins i, i + 1: 

E hl = a, E ia = 1 - c\c % . (1) 

The elementary operators generating the dynamics are 
given by 

* = (°-+ + "m)(*r-^+i). ( 2 ) 

where a i — (of ± icrf)/2 are spin-1/2 raising and low- 
ering operators acting on spin i. In the continuous time 
limit, the operators Cj correspond precisely to Lindblad 
quantum jump operators and generate a dissipative evo- 
lution described by a quantum master equation, cf. Ap- 
pendix]^] The operators act bi-locally on pairs of spins, 
as visualized in Fig. [T|d: Physically, they map any an- 
tisymmetric component in the wave function on a pair 
of sites into the symmetric one, or - in the language 
of hardcore bosons - symmetrically delocalize particles 
over pairs of neighboring sites. Since this process takes 
place on each pair of neighboring sites, eventually only 
the symmetric superposition of spin excitations over the 
whole array persists as the stationary state of the evolu- 
tion: Iteration of the dissipative dynamical map attracts 
the system towards a unique dynamical fixed point, or 
dark state, characterized by p{tg +1 ) = p(ti) = pn, re- 
sulting from the property Cipp = for all i separately. 
More specifically, for m spin excitations initially present 
in the array of N spins, this pure dark state is given by 
the Dicke state 

JV 

\i>) D = \D(m,N))~(j2vt) \D® N 
«=i 

with m collective spin excitations. The derealization of 
the spin excitations over the whole array gives rise to the 
creation of entanglement and quantum mechanical off- 
diagonal long range order, witnessed, e.g., by the single- 
particle correlations (ct+<7~) ^ for \i — j\ — > oo (see 
Appendix [C]) . 

In our interacting lattice spin system, competing uni- 
tary dynamics can be achieved by the stroboscopic re- 
alization of coherent maps p(t£) i— > p(ti + \) = U p{ti)U^ 
with U — exp(— i<j)H) according to the dimensionless spin 
Hamiltonian 

H = ItXtli ® ItXtki = £(! + + of+i)/4. 

i i 

The bi-local terms of the Hamiltonian describe interac- 
tions of spin excitations or hardcore bosons located on 
neighboring sites (see Fig.JlJ,). The competition between 
dissipative and unitary dynamics arises since the dissipa- 
tive dark states pu are not eigenstates of the Hamilto- 
nian, which is diagonal in Fock space and thus leads to 
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Figure 1. Competing dissipative and Hamiltonian dynamical maps in the spin or hardcore boson model, a, 

Dynamical maps, acting on the reduced density matrix p of an open many-body quantum system, can be composed of elementary 
(i) dissipative and (ii) unitary maps. The dissipative maps considered here drive the system into a pure, long-range ordered 
many-body "dark" state po = \iP)d{iP\- The addition of suitable Hamiltonian maps leads to a competition of dissipative and 
coherent dynamics, in such a way that for strong enough coherent interactions, a phase transition into a mixed, disordered 
state results in large systems, b, Analogy of spins and hardcore bosons, in turn equivalent to bosons in the limit of low lattice 
filling (cf. Appendix El: In the considered open-system spin- 1/2 (hardcore boson) lattice model, spin excitations |"|") (||)) can 
be identified with occupied (empty) lattices sites. Quasi-local dissipative maps Dij+i acting on neighboring pairs of spins lead 
to derealization of spin excitations (hardcore bosons) over pairs of lattice sites. Competing coherent dynamics is realized by 
unitary maps corresponding to interactions of spin excitations (hardcore bosons) located on neighboring lattice sites Ujj+i. c, 
Dynamics for an open-boundary chain of N = 10 spins, starting in an initial state containing m = 3 spin excitations (hardcore 
bosons) (see Appendix [5]) . Time evolution is realized by sequentially applying composite maps consisting of (i) dissipative 
and (ii) coherent elementary dynamical maps, as shown in a. For vanishing Hamiltonian interactions, perfect long-range order 
(as measured by the overlap fidelity with the Dicke state \D(m, N))) is built up by the sequential application of dissipative 
maps (white), whereas increasing competing interactions (from yellow to red) lead to a decrease in the dissipatively created 
long-range order. The lower part of the figure shows a zoom into the composite maps 14 to 18, resolving that each of them is 
built up from 9 two-spin elementary dissipative maps, followed by 9 two-spin coherent maps, the latter which can be realized 
by a single, global unitary map in our setting. 



a dephasing of the dissipatively induced off-diagonal or- 
der. The value of the angle <j) S [0, tt/2] determines the 
strength of the competition between the Hamiltonian and 
dissipative dynamics and plays a role analogous to a di- 
mensionless ratio of energy scales in a quantum phase 
transition, or of an energy scale and a dissipative rate 
in the above scenario. Clearly, for small system sizes, 
the sharp transition found in the thermodynamic limit is 
replaced by a smooth crossover as indicated in Fig. [l]:. 



III. EXPERIMENTAL REALIZATION 

The simulation is performed with Ca + ions, confined 
to a string by a macroscopic linear Paul trap (see Ap- 
pendix [C]). Each ion hosts a qubit or spin-1/2, which is 
encoded in the 45i/ 2 (m = -1/2) = |1) = |t) and the 
3L» 5/2 (m = -1/2) = |0) = ||) states. The backbone 
of this digital quantum simulator setup is a universal set 
of high-fidelity operations, which are realized by exactly 
timed laser pulses resonant with the qubit transition (see 
Appendix [C]) . The entangling gate operations [35] act on 
the entire ion string, but the elementary dissipative maps 
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D{ t i + i act on only two of the N system spins (Fig. [2k). 
We achieve operations on subsets of ions via decoupling 
all ions not involved in the elementary map, by shelving 
their population into additional storage states (Fig. ^Ijp) . 
In these electronic states, these ions are effectively "in- 
active" as they do not interact with the globally applied 
laser beams. This novel spectroscopic decoupling tech- 
nique is experimentally simpler than physically moving 
the ions with respect to the laser beam |4"9"] . 

To observe the complex dynamics of the open inter- 
acting spin system, we combine these experimental tech- 
niques (i) to generate long-range phase coherence of spins 
by an iteration of purely dissipative maps, (ii) to com- 
bine these dissipative dynamics with competing coher- 
ent maps, (iii) and finally to implement quantum non- 
demolition (QND) readout and quantum-feedback pro- 
tocols for error detection and stabilization. 

(i) The basis of the composite maps is a single ele- 
mentary dissipative map -D^i+i that is implemented by 
a quantum circuit of coherent gate operations and ad- 
dressed optical pumping (see Fig. 2fc), acting on the 
two currently active ions i, i + 1 and an ancillary qubit, 
which is used to engineer the coupling to the environ- 
ment [T31 Q31 [55] . The circuit decomposition of the three- 
qubit unitary underlying a single elementary dissipative 
map, is obtained from an optimal control algorithm (see 
Appendix C2|, resulting in a sequence of 17 operations 
containing 4 entangling gates. We have characterized 
a single elementary dissipative map by quantum pro- 
cess tomography on the two system qubits leading to a 
mean state fidelity of 68(1)%, which approximately corre- 
sponds to an average fidelity of over 98% per gate opera- 
tion (see Appendix |C]). Due to the considerable complex- 
ity of the gate sequence, errors occur in different bases 
and are expected to average out and give rise to depo- 
larizing noise without any preferred direction. There- 
fore, the actual implemented dynamics can be modeled 
as a combination of the ideal, tailored dissipative map 
and the depolarizing noise channel. Detailed numeri- 
cal simulations show that the error is mainly caused by 
laser frequency and magnetic field fluctuations (see Ap- 
pendix IE 2h . 
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We repeatedly apply such elementary maps D.^+i, in- 
terspersed with decoupling pulses to coherently transfer 
ions in and out from the storage states, to build up the 
composite dissipative map in a modular way, Fig. [2^,. Us- 
ing 3+1 ions, we studied pumping towards Dicke states in 
a three-spin chain with open boundary conditions, where 
we applied up to three simulation timesteps, each con- 
sisting of the two elementary maps D\ 2 and D 2 ^. The 
results displayed in Fig. [3^, show a clear experimental sig- 
nature of dissipatively induced derealization of the spin 
excitations during the application of the first two elemen- 
tary dissipative maps, before experimental imperfections 
become dominant for longer sequences. 

(ii) To investigate the competition between dissipative 
and competing Hamiltonian dynamics, we added elemen- 
tary unitary maps according to nearest-neighbor spin- 




Figure 2. Experimental procedure to implement 
open-system dynamical maps, a, Schematic overview of 
the experimental implementation of a composite dynamical 
map consisting of (i) multiple elementary dissipative maps, 
(ii) coherent competition, and (iii) error detection and cor- 
rection. Decoupled ions are represented as gray bullets and 
decoupling (re-coupling) operations as gray (blue) squares, b, 
Scheme for decoupling ions from the interaction with the ma- 
nipulating light fields: (i) shelve population from 451/2 (m = 
-1/2) = |t) to 3D 5/2 (m = -5/2), (ii) transfer the popu- 
lation from 3D 5/2 (m = -1/2) = ||) to 4S 1/2 (m = +1/2), 
and subsequently to (iii) 3D 5 / 2 (m = —3/2). c, A single dis- 
sipative element is realized using two system spins and one 
ancilla qubit (|0) = |4},|1) = |t}) by (i) mapping the informa- 
tion whether the system is in the symmetric or antisymmet- 
ric subspace onto the logical states |1) or |0) of the ancilla, 
respectively; (ii) mapping the antisymmetric onto the sym- 
metric state using a controlled phase flip conditioned by the 
state of the ancilla qubit; and finally (iii) reinitialization of 
the ancilla qubit via optical pumping using the ^ 2 P\/2 state 
(see Appendix C4l. d, Schematic view of the competing in- 
teraction consisting of quasi-local unitary maps Ujj+i. 



spin interactions to the composite dynamical maps (see 
Appendix [D]). Note that due to the commutativity of 
the two-spin elementary unitary maps, the composite, 
globally acting unitary map can be realized by a sin- 
gle unitary sequence (Fig. [2ji) . The results displayed in 
Fig. [4] for experiments with 3+1 and 4+1 ions show a 
clear fingerprint of incompatible Hamiltonian dynamics, 
which competes with the dissipative maps driving the 
spin chains towards the Dicke states. Further measure- 
ments with varying excitation number and competition 
strength are discussed (see Appendix [H| . 

(iii) Under the ideal dynamical maps, the total spin 
excitation (or hardcore boson) number m is a conserved 
quantity, however the presence of depolarizing noise as 
the dominant error source results in population leakage 
out of the initial excitation number subspace. To re- 
duce this detrimental effect and enable the implementa- 
tion of longer sequences of dynamical maps, we devel- 
oped and implemented two counter-strategies (see Fig. [5] 
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for details). In a first approach we applied a quantum 
non-demolition (QND) measurement of the spin excita- 
tion number at the end of the sequence of dynamical 
maps, which allowed us to detect and discard experimen- 
tal runs with a final erroneous excitation number and 
thereby improve the overall simulation accuracy. This 
global measurement is QND in the sense that only infor- 
mation about the total number of excitations, but not 
on their individual spatial locations along the chain is 
acquired; thus the simulation subspace is not disturbed 
(see Appendix [f]). The results shown in Fig.[3]and Fig. [4] 
confirm that the errors in the spin excitation number arc 
strongly suppressed and a reasonable overlap with the 
ideal evolution can be maintained for more simulation 
time steps. 

Complementary to this post-selective method, we in- 
troduced a more powerful, active QND feedback scheme, 
which bears similarities to quantum feedback protocols as 
realized with photons in a cavity [50] ■ The key idea is to 
actively stabilize the spin system during the sequence of 
dynamical maps in a subspace with a particular spin exci- 
tation (or hardcore boson) number (see Appendix [Gil. In 
order to be able to perform this stabilization with a sin- 
gle ancilla qubit, we break the stabilization process into 
two independent parts, where the first part removes one 
excitation if there are too many excitations in the sys- 
tem, and the second part adds one excitation if needed. 
Similarly to the post-selective technique presented above, 



first the information whether there are too many (few) 
excitations in the system is coherently mapped onto the 
ancilla qubit. Depending on the state of the ancilla qubit, 
a single excitation is removed from (injected into) the 
system by a quantum feedback protocol. This extraction 
(injection) is in general an ambiguous process, as the ex- 
citation can be removed (injected) on multiple sites. We 
use a scheme that tries to perform the removal (injec- 
tion) subsequently on each site and stops once it was 
successful. Using only a single ancilla qubit, this pro- 
cess cannot be performed efficiently as a unitary process, 
therefore we developed a technique making use of the re- 
setting and decoupling techniques described above (see 
Appendix [G|. 

We demonstrate the excitation removal for a chain of 
3+1 spins, initially prepared in an equal superposition of 
all basis states, as shown in Fig. [3]:. At the current level 
of experimental accuracy, the implementation of this sta- 
bilization scheme cannot improve the performance when 
used in the full simulation sequence (see Appendix |G 4[ ) . 
We emphasize, however, that our approach relies only 
on a single ancillary qubit, regardless of the system size. 
More generally, such customized error detection and re- 
duction strategies will incur a substantially reduced re- 
source overhead as compared with full-fledged quantum 
error correction protocols. For the future, we therefore 
envision that methods as those demonstrated here will 
be a promising practical avenue to considerably extend 
the runtime of quantum simulators. 
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Figure 3. Experimental results of dissipatively in- 
duced delocalization through composite dynamical 
maps with 3+1 ions. The results from an ideal model 
are shown in light-blue bars whereas those from a model in- 
cluding depolarization noise are indicated by dark-grey bars. 
Blue rectangles indicate the experimentally observed dynam- 
ics without any correction scheme whereas red diamonds in- 
clude a post-selective error detection scheme (error bars, la), 
a, Dissipative pumping into a three-spin Dicke state: Starting 
in an initial product state with two localized spin excitations 
| tit) i the application of the first two elementary dissipative 
maps leads to an increase in the delocalization of the two exci- 
tations over the spin chain, which is reflected by an increasing 
state overlap fidelity with the three-spin Dicke state \D(2, 3)}. 
However, after applying a second and a third composite dissi- 
pative map, a decrease in the state overlap fidelity sets in and 
becomes dominant for long sequences of dynamical maps, b, 
The presence of depolarizing noise results in population leak- 
age out of the initial subspace with m = 2 spin excitations. 
This effect is evident in the decay of the probability of find- 
ing the three-spin system in the m = 2 excitations subspace 
as a function of the number of applied elementary dissipative 
maps. A single composite dissipative map is indicated by a 
yellow bar. 
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Figure 4. Experimental results for competing dissipa- 
tive and coherent dynamics with 3+1 and 4+1 ions. 

As in Fig. [3] the results from an ideal model are shown in 
light-blue bars whereas those from a model including depo- 
larization noise are indicated by dark-grey bars. Blue rectan- 
gles indicate the experimentally observed dynamics without 
any correction scheme whereas red diamonds include a post- 
selective error detection scheme (error bars, la). The appli- 
cation of dissipative (coherent) maps is indicated by yellow 
(red) bars. Competing dissipative and coherent dynamics for 
m — 2 excitations in chains of a, N = 3 and b, N = 4 spins: 
the spin chains are first driven towards the Dicke-type dark 
state by the two and three elementary dissipative maps for a 
system of 3 and 4 spins. The subsequent application of the 
non-compatible unitary dynamical maps leads to a strong de- 
crease of the overlap with the respective Dicke states, before 
subsequent elementary dissipative maps again start to pump 
the system back towards the Dicke states. 
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Figure 5. Experimental error detection and reduction techniques, a, Numerical simulation of the time evolution con- 
sisting of composite dynamical maps building up long-range order quantified by overlap fidelity with the Dicke state \D(m, N)), 
starting in an initial state containing m = 3 spin excitations on an open-boundary chain of N — 10 spins. The simulation 
was performed for (i) ideal operations without experimental imperfections, (ii) a model including experimental errors, (iii) a 
model including errors but with the stabilization scheme. The lower part depicts the expected excitation distribution of the 
time evolution demonstrating that the stabilization keeps the system in the correct excitation subspace. b, Schematic idea 
and quantum circuit for a post-selective QND measurement of the total spin excitation number: First the information whether 
the N system spins are (not) in the subspace with m spin excitations or not is mapped coherently onto the logical state |1) 
(|0)) of an ancillary qubit. A subsequent projective measurement of the latter indicates the presence of an erroneous excitation 
number, in which case the experimental run is discarded, c, Active QND feedback scheme to stabilize the spin system in a 
desired excitation number subspace by actively extracting (injecting) spin excitations in case errors in previous dynamical maps 
have led to a larger (smaller) excitation number than present initially. The information whether more excitations than the 
desired value are present in the system (or not) is coherently mapped onto the state |1) (|0)) of the ancilla qubit. A non-unitary 
controlled-operation, only active for the ancilla qubit in |1), then extracts in a minimally destructive way one spin excitation 
from the system and automatically stops once this is achieved, d, Experimental demonstration of the stabilization protocol 
for mo = 1 using 3+1 ions by applying the excitation removal to an initial state consisting of an equal superposition of all 
computational states, \ipo) = 1 / a/8( i) + | t)) 83 - Here, a single spin excitation should be removed if two or more excitations are 
present in the system (blue bars). Thus, pumping from the m — 3 into the m — 2 and from the m = 2 into the m = 1 subspace 
is expected, whereas population initially present in the m — and m = 1 subspaces should be left untouched. The correct 
populations after the application of the excitation removal protocol are confirmed by the measured data (blue rectangles) which 
are close to the ideally expected behavior (colored bars). Data error bars are smaller than the markers. 
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These appendices provide background material on de- 
tails of the theory and the experimental implementation 
of the quantum simulation of open system dynamical 
maps with trapped ions discussed in the main text. In ad- 



dition, it covers further experimental simulation results 
which are not discussed in the main text. Appendix [A] 
discusses how the continuous time evolution of a Lind- 
blad master equation is recovered in the limit of infinites- 
imal time steps of Kraus maps. Appendix |TT| shows how 
the dissipative spin operators are constructed in analogy 
to the case of bosons. Appendix [C] provides details on 
the experimental system and the available set of tools 
for open-system quantum simulation. In Appendix [T5] 
we present in more detail the general concept of the en- 
gineering of elementary dissipative and coherent Kraus 
maps, and also provide the specific sequences of opera- 
tions that have been used in the experimental implemen- 
tation. In Appendix[E]we discuss the quantitative exper- 
imental characterization of elementary dissipative Kraus 
maps. Based on these measurements, we describe in more 
detail a theoretical noise model, which captures the dom- 
inant experimental imperfections, and present analytical 
expressions and information on the numerical simulations 
for the dynamics of the open many-body spin system un- 
der dynamical maps. Appendices [F] and [G] contain back- 
ground information on the error detection and reduction 
schemes, as well as experimental details on their imple- 
mentation. Finally, in Appendix [H] we discuss additional 
experimental simulation results, including the realization 
of competing dissipative and coherent dynamical maps 
for varying Hamiltonian interaction strengths and initial 
states, and the combined implementation of dissipative 
dynamical maps with the error reduction protocol. 

Table |TJ provides an overview of the algorithmic build- 
ing blocks required for the implementation of elementary 
dissipative Kraus maps, the realization of Hamiltonian 
dynamical maps, spectroscopic decoupling of ions, and 
elements for the detection and reduction of errors in our 
ion-trap quantum simulator. The table describes the ex- 
perimental resources needed for the implementation of 
these elementary sub-routines as well as for the combina- 
tion of these building blocks in more complex composite 
dynamical maps. The details of the underlying Kraus 
map engineering as well as the detailed quantum circuits 
used in the experiment are presented in the correspond- 
ing appendices below. 



Appendix A: Dynamical maps vs. quantum master 
equation 

The dissipative Kraus maps specified in equation ([TJ 
are obtained as a special case of the operators En — 
sin0Ci, E t _ 2 = 1 + (cos6> - l)c\ci, for 9 = n/2. This 
limit corresponds to a deterministic action of the Kraus 
map (see Appendix II), in this case generating truly 
stroboscopic dynamics. Instead, in the limit 9 — > 0, 
we approximate En ~ 0a, En « 1 — §0 2 cjci. In 
this limit, the sequence of dynamical maps reduces to 
the continuous time evolution described by a quantum 
master equation entirely generated by a dissipative Liou- 
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5 


4 
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C3 




Composite dissipative map (3 spins) 


12 


26 


24 


2 


64 




Composite diss, and coh. dyn. map (3 spins) 


14 


29 


26 


2 


71 




Composite diss, and coh. map + QND (3 spins) 


18 


40 


27 


3 


88 




Composite dissipative map (4 spins) 


18 


33 


35 


3 


89 




Composite diss, and coh. dyn. map (4 spins) 


21 


44 


44 


3 


112 





TABLE I. Summary of the required resources for the elementary and composite dynamical maps and additional tools used 
in the quantum simulation. The required operations for the composite maps do not strictly match the sum of the required 
elementary operations, since in the implementation of composite maps synergy effects in the resources for the spectroscopic 
decoupling operations can be exploited. 



ville operator, £[p] = \CiP c t ~ H c l c i>P}) P3' Sim- 
ilarly, the Hamiltonian Kraus map can be expanded, 
cxp(— i<fiH) pa 1 — i<pH for <j> — > 0. The dynamics in 
this continuum limit is then described by the quantum 
master equation dtp = — i[UH, p] + k£[/>], with the di- 
mensioned energy scale U (Udt = <j)) and dissipative 
rate k (ndt = 2 ). Here, dt is the physical time required 
for the implementation of one Kraus map in the digital 
simulation. The dimensionless ratio describing the com- 
petition g = U/k — 4>/9 2 remains well-defined in this 
limit. We finally note that a temporal coarse graining 
implemented by averaging over a sequence of elementary 
maps with even larger discrete mapping steps gives rise 
to a quasi-continuous evolution, as shown numerically in 
Fig. [Tj:. In a large system, the quasilocal operations can 
be coarse grained also over space. The result is an effec- 
tive quasi-continuous master equation dynamics for the 
density operator. 

Appendix B: Atomic boson vs. ionic spin model 

The dissipative spin operators Ci of equation ([2| are 
constructed in complete analogy to the case of bosons, 
which has been proposed theoretically in [29) : Formally, 
and as further detailed in Appendix [Dj they obtain by 
replacing the spin raising (lowering) operators af (a^ ) by 
bosonic creation (annihilation) operators a\(a.i) °f atoms 
confined to an optical lattice. In that case, the dark state 
is given by m symmetrically delocalized particles on N 

lattice sites, i.e. \^) D = (ml)- 1 / 2 fe^al)"* \®)® N ~ 
the Dicke dark states of the spin model are replaced by a 
fixed number Bose-Einstein condensate (the bosonic vac- 
uum state is defined by cii\0)® N = for all i). Using the 



Holstein-Primakoff representation of spin 1/2 operators 
in terms of bosons, af = a\\/l — hi (hi = a]<Xj), it is 
seen that the dissipative spin operators reduce to their 
bosonic counterpart in the limit of small average occupa- 
tion h = (hi) <C 1, where the square root can be safely 
replaced by one. 

Appendix C: Experimental system and techniques 

In this appendix we will describe the available opera- 
tions in our universal ion-trap quantum simulator. 

1. Coherent gates 

The qubit is encoded in the 4«S , 1 / 2 (fTij = —1/2) and 
3D 5 / 2 (mj = —1/2) electronic states of the 40 Ca + ion 
and is manipulated by precisely timed light pulses on res- 
onance or near resonant with the optical transition. The 
laser light can be applied from two different directions as 
depicted in figure |A1[ where one beam illuminates the 
entire ion string homogeneously and the second beam is 
able to address each ion individually |A1| . 

The set of coherent gates consists of collective single- 
qubit rotations, addressed single-qubit gates and collec- 
tive entangling gates. Collective single-qubit rotations 
are implemented by the globally applied laser beam, res- 
onant with the qubit transition, realizing the unitary 

R(e,<f>) = exp(- l ^S 4 ), (CI) 

with Stj, = Yli=o a t a collective spin operator, and 
af — <Ji cos(0 • 7r) + af sin(0 • tt) 
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Figure Al. Illustration of the geometry of the laser-beams 
used for qubit manipulation. A global beam illuminating the 
ion string homogeneously is used to implement collective ro- 
tations R(6, (j>) and multi-ion M0lmer-S0rensen-type entan- 
gling gates MS (8, (/>), whereas an addressed beam enables the 
realization of single-qubit rotations Sz(9,j). An addressed 
optical pumping technique allows one to incoherently reini- 
tialize the state of an ancillary qubit, which plays the role of 
a tailored environment. See text for more details. 



being a linear combination of the single-qubit Pauli ma- 
trices af and of, acting on qubit i. The rotation angle 
9 is determined by 9 = Qt/tt, which can be controlled 
by the Rabi frequency f2 and the pulse duration r. In 
this notation a complete 7r-flop inverting the electronic 
population of the logical states corresponds to a rotation 
angle 9 = 1. On the other hand, tuning the phase <j> of 
the global laser beam allows one to control the rotation 
axis lying in the x-y-plane of the Bloch sphere, around 
which each of the qubits is rotated. Addressed single- 
qubit operations 



S Z (V, i) = cxp ( — —a\ 



(C2) 



around the z-axis are realized by shining in laser light 
near resonant with the qubit transition that induces an 
intensity-dependent AC-Stark shift Aac- Again, the 
rotation angle 9 is determined by the pulse length r, 
9 = Aact/k, and a 7r-pulse corresponds to 6 = 1. Fi- 
nally, collective entangling operations are implemented 
by a bi-chromatic, globally applied laser field, which ef- 
fectively realizes two-body M0lmer-S0rensen (MS) type 
interactions, 



MS(9, 0) = exp 



i9n 



exp 




between all pairs i and j of the ion chain — 
0, 1, ... N) gHl El]. Again, the angle 4> allows one to 
control whether erferj (for cf) = 0), of crj' (for <p = 1/2) 

or interactions o'faj corresponding to any other axis 
in the x-y-plane are realized. In this notation, the an- 
gle 9 = 1/2 corresponds to a "fully-entangling" MS gate, 



i.e. a unitary which maps the computational basis states 
of + 1 ions onto multi-particle entangled states, which 
are (up to local rotations) equivalent t o N + 1-qubit GHZ 
states. Altogether, operations in Eq. (CI) to (C3) form 



a universal set of gates, enabling the implementation of 
arbitrary unitaries on any subset of ions |A3| . 



2. Numerical optimization of gate sequences 

Any unitary operation required for the implementa- 
tion of dissipative and coherent maps, as well the error 
detection and reduction protocols, needs to be decom- 
posed into a sequence of available operations. As dis- 
cussed in more detail below, such decompositions can be 
constructed systematically |A4j . However, as such de- 
compositions are in many cases not optimal in terms of 
the number of required gate operations, it is convenient 
to resort to a numerical optimal control algorithm |A3j 
to search for optimized sequence decompositions involv- 
ing less gates. Whereas the numerical optimization algo- 
rithm becomes inefficient for general unitary operations 
acting on a large number of qubits, it is well-suited for 
the optimization of unitaries which act only on a small 
subset of ions (such as 2+1 ions in the implementation 
of an elementary dissipative Kraus map), independently 
of the total system size. 

Numerically optimized pulse sequences may include 
global AC-Stark pulses and MS gates with negative ro- 
tation angles 9 < 0, which are not directly contained in 
the available gate set discussed in Sec. |C 1| However, 
as any collective rotation around the z-axis of the Bloch 
sphere can be interpreted as a re-definition of the x- and 
y-axes, a global AC-Stark pulse can be omitted if the 
phases <f> of the following resonant operations are prop- 
erly adjusted. Regarding MS gates with 9 < 0, these can 
be implemented by MS gates with positive rotation an- 
gles, as MS(—9, <p) = MS(1 — 9, 4>) up to local rotations 
(see Eq. (9) in Ref. 



3. Spectroscopic decoupling of ions 

Despite the globally applied beams for the collective 
rotations and MS gates, operations on subsets of ions 
can be realized by spectroscopically decoupling ions not 
involved in the realization of a certain Kraus map from 
the dynamics. This is realized by coherently transferring 
ions to electronic states, where they do not couple to the 
globally applied light fields to a very good approximation. 
This decoupling technique enables the use of optimized 
sequences for the realization of Kraus maps on a small 
number of sites, independently of which subset of ions is 
currently involved in the map, and independently of the 
system size, i.e. the total number of sites. 

The full decoupling sequence consists of the following 
parts as also outlined in Fig. 2 in the main-text: (i) 
First, the population is transferred from the qubit state 
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Transition 


Pulse 


\S) -> \D') 
\S) \D') 
\S) -> \D) 
\S) -> \D') 


i?(0.5,0) 
S z (l,i) 

s,(i,j) 

#(0.5, 1) 


\D)^\S>) 
\D) ->■ |S") 
|D) ->■ S") 
|D)^|S') 


i?(0.5,0) 
5,(1, i) 
3.(1,3) 

R(0.5, 1) 


|S") ->■ |D") 


R(i,o) 



TABLE II. Pulse sequence for spectroscopic decoupling of ions 
i and j by coherently transferring their quantum information 
from the qubit states |S) and \D) to the states \D") and 
\D"). The gates listed in the three blocks of the table realize 
the pulses (i) - (iii) shown in Fig. 2b of the main text. 



\S) = 4S 1/2 (m / = -1/2) = |1) = | t) to the \D>) = 
3D 5 / 2 (fnf — —5/2) state, (ii) Then, the population from 
the other qubit state \D) = 3D 5/2 (m f = -1/2) = |0) = 
| I) is transferred via |S") = 4S , 1 / 2 ("i/ = +1/2) to the 
\D") = 3D 5 /2(mf = —3/2) state. The required pulses for 
decoupling two ions are shown in Table [IT] Bringing the 
population from the decoupled states back to the origi- 
nal qubit states is realized by implementing the described 
sequence in reverse order. 

The decoupling technique introduces additional errors 
that are not included in the theoretical error model. A 
rigorous treatment of these errors is cumbersome since 
it cannot be modeled in a qubit system anymore, but 
involves the full electronic substructure of the ion. How- 
ever, the effect of the decoupling process in the compu- 
tational basis can be characterized by quantum process 
tomography. We found a process fidelity of 94(2)% for 
decoupling a qubit and bringing it back to the original 
states. Next, we proved that the decoupled qubit is in- 
deed to a high degree not affected by the manipulation 
pulses. We checked this in a system of 3+1 ions, where we 
first decoupled a single qubit, then applied the pulses as 
required for a single elementary dissipative Kraus map 
on the remaining two system ions and the ancilla ion, 
and finally transferred the decoupled qubit back to the 
original qubit state. Due to residual far off-resonant cou- 
pling to transitions coupling different Zeeman sublevels, 
the pulses implementing the dissipative Kraus map in- 
duce a deterministic AC-Stark shift on the decoupled 
ion. This Stark shift is measured with a Ramsey-type 
experiment, and its compensation is performed with the 
final two pulses in the sequence as shown in Table [H| 
Quantum process tomography on the decoupled qubit, 
where the systematic Stark shift has been compensated, 
yields a process fidelity of 93(2)%. Thus we can conclude 
that the pulses corresponding to the Kraus map imple- 
mentation do not affect the decoupled qubit significantly, 
and that the dominant errors result from laser intensity 
fluctuations in the decoupling pulses themselves. 



Transition 


Pulse 


\S') -> \D) 
\S') -> \D) 
\S') -> \D) 


.R(0.5,0) 

Sz(l,j) 

#(0.5,1) 


\S') -> \P) 


a~ repump 



TABLE III. Pulse sequence for the individual reset of qubit j. 



4. Incoherent reinitialization of individual ions 

The implementation of an elementary dissipative 
Kraus map is completed by an incoherent reset of the an- 
cillary ion to its initial state \S), see step (iii) in Fig. 2c 
in the main text. This reinitialization is realized by an 
optical pumping technique: First, an addressed pulse is 
applied to the ancillary ion to transfer the population 
present in the \D) state to the \S') state. Then cr~ po- 
larized light is applied to the entire ion string performing 
optical pumping from \S') towards \S) via the short-lived 
4 2 -Pl/2 state. This procedure does not affect the infor- 
mation in the system ions encoded in the original qubit 
states 15*) and \D), as the light couples only to the \S' 
level The required operations are shown Table III 



Appendix D: Engineerging of dissipative and 
Hamiltonian dynamical maps 

In this appendix we provide details on the engineering 
and the specifics of the circuit-based experimental im- 
plementation of elementary dissipative and Hamiltonian 
dynamical maps. 



1. Action of the dissipative Kraus maps 

The elementary dissipative Kraus maps 

with 

and 



Ei,i = Ci 



Ei,2 - 1 



c]ci 



(Dl) 



(D2) 



are generated by the operators Cj (cf. Eq. (2) in the main 
text), as given by 



(D3) 



These operators act bi-locally, i.e. involve two neighbor- 
ing spins i and i + 1, whereas the other spins are specta- 
tors. It is instructive to examine their action on the basis 
states of the local Hilbert space of the two qubits i and 
i + 1, which is spanned by the singlet and triplet states 
of the total spin Sf i+1 of the two spin-1/2 particles: 

5? )<+ i^(|01)-|10)) = 
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Figure A2. Schematics of the action of the two-spin operator 
(D3l in the local Hilbert space of the two spins i and i + 1, 



which is spanned by the singlet state (left sector of the Hilbert 
space) and the triplet states (right sector). The operator a 
converts the singlet state +g (|01) — 1 10}) into the triplet ms = 

state ^75 (1 01) + 1 1 0) ) . All triplet states are dark states and 
left invariant. 



and 



S* i+1 (|00), ^(|01) + 110)), |11)) 
= 2(|00), ~(|01) + |10)), |11)). 



Here, S i)i+ i = S, t + S, i+1 with S l 
and 



For simplicity of the notation, we skip the spin indices i 
and i + 1 for the states and use the short-hand notation 



1 00) = |0)i ® |0) i+ i, etc. As illustrated in Fig. [A2j the 
operators Cj induce pumping from the singlet into the 
triplet ms = state, 

Cj -L(|oi)-|io» = -^(|oi) + |io)) 1 

whereas all triplet states are dark states, 

1 00) = Ci-^(|01) + |X0>) = a\ll) = 0. 

As shown in Fig. lb of the main text, under this dis- 
sipative dynamics a single spin excitation (or hardcore 
boson) is symmetrically delocalized over the two sites, 
whereas the states of two (|11)) or zero ( 1 00) , "vacuum") 
spin excitations (or hardcore bosons) are left unchanged. 



2. Circuit-based implementation of elementary 
dissipative Kraus maps 

Open-system dynami cs a ccording to the elementary 
dissipative Kraus maps ( Dl ) can be realized in a "digital" 



way, by using quantum simulation tools for open systems, 
which have been previously developed and demonstrated 
experimentally in the context of dissipative preparation 



controlled flip of 
spin i (ii) 



mapping (i) 




>0»- 
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• — 
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{exp( 


±i<t>a z )} 


i+1 • — 









inverse 
mapping (iii) 



deterministic flip 
of spin i (ii) 



mapping (i) 

reset of the 
ancilla (iv) 

-CJB-I1) |i>o 



reset of the 
t y ancilla (iii) 




Figure A3. Quantum circuits for the realization of elemen- 
tary dissipative Kraus maps (Dl I on spins i and a, First, 
the information whether the two system spin-1/2 particles i 
and i+1 are in the singlet or triplet subspace of the two-spin 
Hilbert space (see Fig. |A2[ ) is coherently mapped onto the 
logical states |0) and |1) of an ancilla qubit which is initially 
in |1). (ii) Next, a two-qubit gate is applied, which realizes an 
effective spin flip on one of the system spins, and thereby con- 
ditionally on the state of the ancilla converts the singlet into 
the triplet m s = state: ^j(|01) - |10)) -> ^j(|01) + |10». 
After inverting the mapping (iii), the ancilla ion is optically 
pumped backed to its initial state |1) (iv). This last steps 
renders the dynamics irreversible and provides the dissipative 
ingredient to extract entropy from the two system spins. The 
probability of the conversion from the singlet to the triplet 
state is controlled by the angle <j) appearing in the two-qubit 
gate, and it is given by p = sin 2 <f>. In the limit of <j> <C 7r/2, 
the conversion from the singlet into the triplet state takes 
place only with a small probability, and the general dissipa- 
tive Kraus map describing the pumping process reduces to 
the master equation (D8l. b, In this work we are interested 



in deterministic pumping from the singlet into the triplet sub- 
space, i.e. in the case of probability p = 1 for (f> — n/2. In 
this situation, the circuit simplifies as the inverse mapping to 
partially disentangle the system spins from the ancilla is not 
required, and the circuit simplifies to the three-step process, 
which is also shown in Fig. 2c of the main text. 



of Bell and multi-qubit stabilizer states [25] . The key idea 
in engineering the two-spin dissipative dynamics accord- 
ing to Eq. (Dl I is to combine the experimentally available 
gates (see Sec. C 1 below for details) with optical pump- 
ing on an additional ancillary qubit (see Sec. C4), which 
plays the role of a tailored environment. 



General engineering strategy - The observation that 
the singlet state is dissipatively converted into the ms — 

triplet state suggests the following gate-based imple- 
mentation via a four-step process (shown in Figure A3 1) , 
which involves a circuit of unitaries applied to the qubits 
i, i + 1 (steps (i) to (iii)), followed by the incoherent reset 
of the ancilla qubit (step (iv)): 

(i) First, a unitary M, acting on the two system spins 

1 and i + 1 and the ancillary qubit, coherently maps the 
binary information whether the two system spins are in 
the singlet or triplet subspace onto the two logical states 
|0) and |1) of the ancillary qubit. This is achieved by the 
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unitary 



M{6) = cxp (- 1 (Si i+1 - 2) ® (D4) 

= exp (i0P M+i ® do ) 

= cos(9P iji+ i) ® lo - isin(0P ijJ+ i) ® (Tq 

for which we choose 9 — tt/2 so that the unitary reduces 
to 



M (tt/2) = (1 - P 



M+JJ 



1 - i Pj,i+i <8> Co 



(D5) 



The unitary M is constructed in a way that the state of 
the ancilla qubit is rotated conditional on the state of 
the two system spins, and that the angle of this rotation 
depends on the operator (Sf i+1 — 2) which acts on the 

two system spins i and i + Here, Pj,j+i = c\ci is 
the projector onto the singlet subspace, (1 — Pi.i+i) its 
orthogonal complement, and we have used Sf j+1 = 2(1 — 
P,i+i)- Under the unitary M(tt/2), the ancilla qubit is 
rotated from its initial state |1) into |0) if and only if the 
two systems spins are in the singlet state. 

(ii) Next, the transfer of the system qubits from the 
singlet to the triplet subspace is enabled by a two-qubit 
gate, which acts on the ancilla qubit and the system qubit 
i. It reads 



C(4>) =exp 

=|1)<1|o®1< + 



4 a 



55 exp(±j$of ). 



(D6) 



This controlled-operation applies a z to the i-th system 
qubit, and thus with probability p = sin 2 4> converts 
the singlet into the triplet state: 0f4g(|Ol) — 1 10) ) = 

4g (|01) + |10)). Here, the ancilla qubit acts as a quantum 

controller [141 116j , which is not observed but controls the 
coherent feedback which is applied to the system qubits. 
From Eq. (D6) it is clear that if the ancilla is in state 
|1), corresponding to the case of the two system spins 
residing in the triplet subspace, no quantum feedback is 
applied and the state of the system spins remains un- 
changed. Due to the previous mapping M, this assures 
that the triplet subspace is left invariant. 

(iii) The initial mapping (i) is inverted by applying 
the inverse unitary AP(7r/2). For the system being in 
the target triplet space from the outset, this re-installs 
the conditions before step (i) . For the special case of de- 
terministic pumping (p — 1), as experimentally realized 
in this work, the inverse unitary can be omitted (see also 
Figure A3 a). 



(iv) After these unitary steps, the ancillary qubit re- 
mains in general entangled with the two system spins i 
and i + 1. It is finally reset to its initial state |1) by 
optical pumping: this is the physical dissipative mecha- 
nism, which renders the dynamics irreversible and allows 
to extract entropy from the system spins. This refreshes 
the ancilla qubit and prepares it in the known pure state 
|1), so that it can be used for the implementation of sub- 
sequent dissipative Kraus maps. 



For an initially uncorrelated state of the ancilla qubit 
and the system spins, |l)(l|o €5 p, the resulting dynamics 
for the system spins (described by the reduced density 
matrix p of the system spins) is obtained by applying 
the combined unitary U = M^(tt/2) C{4>) M(tt/2) and 
tracing over the ancilla qubit's degrees of freedom, 



p^tr {U(\l)(l\ ®p)U^}. 



(D7) 



Straightforward algebra yields the Kraus map p h-» 
£fc=i 2 Ei,kpEj u with the operation elements 



E it i = sin^Ci 



and 



Pi. 2 = 1 + (cos 0-1) c\ci 



One the one hand, for <f> — tt/2 the operation elements re- 



duce to Eq. (D2| and correspond to deterministic pump- 
ing (p = 1) from the singlet into the triplet states, which 
is the scenario realized in this work. On the other hand, 
for (j) <C 7r/2 one can expand the operation elements of 
the Kraus map and recovers the master equation limit 



P ^ P + 4? (cipct - i {cJ Ci) p} j + O(0 4 ). 



(D8) 



Additional remarks - We note that our implementa- 
tion corresponds to an open-loop control scenario, where 
the ancillary qubit remains unobserved during the sim- 
ulation. However, we remark that it is possible to mea- 
sure the state of the ancillary qubit in the computational 
basis by an addressed fluorescence measurement, before 
it is reset to its initial state by optical pumping [25] . 
Such measurement reveals whether under the application 
a particular dissipative Kraus map the two system spins 
have undergone a collective two-spin quantum jump from 
the singlet into the triplet state or not. This informa- 
tion yields "in-situ" information about the many-body 
system's dynamics along a particular trajectory. The 
temporal statistics of quantum jumps in open, driven 
many-body quantum spin systems contains valuable in- 
formation about dynamical phase transitions |A5[ IA6| . 
In the present work we do not further explore this pos- 
sibility. We remark that the fluorescence measurement 
of the ancillary ion, which is associated to the scattering 
of many photons and non- negligible heating of the vi- 
brational modes, can be combined with subsequent laser 
re-cooling of the relevant vibrational modes via the an- 
cillary ion, as demonstrated by |A7j . This allows one 
to perform such quantum jump measurements and to 
afterwards re-initialize the external degrees of freedom, 
thereby enabling that in the same experimental run fur- 
ther gate operations required for the implementation of 
subsequent Kraus maps can be still applied with high 
fidelity. 

Specific experimental implementation - For the imple- 
mentation of the elementary dissipative Kraus maps in 
the ion-trap simulator, the unitary operations M(tt/2) 
and C(tt/2) (see Eqs. §D4§ and jp6| ) have to be de- 
composed into a quantum circuit of available gates: The 
mapping M can be written as a product of three unitaries 
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involving 3-body spin interaction terms, 
M(?r/2) = exp 



2 - 2' 



(D9) 



exp 



o ex P 



* 7r „x„x„x 
~~8 



exp 



i+l 



'i+l 



As discussed in Ref. |A4j each of the three-body uni- 
taries can be realized by two "fully-entangling" MS gates 
in combination with single-qubit rotations, such that M 
could be implemented by a quantum circuit involving six 
MS gates and a number of single-qubit gates to rotate 
the system spins between the x, y and z-bases between 
the different unitaries. 

Similarly, the two-qubit controlled operation C(tt/2) 
(see Eq. (D6|) is up to local unitaries equivalent to one 
"fully-entangling" two-ion MS gate |A4j . 

C(tt/2) ~ ex P H(^/4Kaf), 

acting on the ancilla qubit and the system spin i. In 
the experiment, we did not use such a systematically 
constructed circuit decomposition, but instead resorted 
to the numerical optimization algorithm described in 



Sec. C 2 to further reduce the complexity of the quan- 
tum circuit: The gates of the experimentally employed 
sequence decomposition for the implementation of the 
unitaries of one elementary dissipative Kraus map are 
listed in Table HVI 



Number 


Pulse 


Number 


Pulse 


1 


5,(1.5,0) 


11 


MS(0.5,1.5) 


2 


#(1.5,1.0) 


12 


5.(1.75,2) 


3 


M5(0.25, 1.0) 


13 


M5(0.5,2.25) 


4 


5 Z (1.0,1) 


14 


#(0.5,1.75) 


5 


M5(0.875, 1.0) 


15 


#(0.5,2.25) 


6 


5,(1.0,2) 


16 


M5(0.25,2.25) 


7 


M5(0.125, 1.0) 


17 


5* (1.5, 2) 


8 


5,(1.5,2) 


18 


5,(1.0,1) 


9 


#(0.5,0.5) 


19 


#(0.5,2.25) 


10 


5,(0.5,2) 







TABLE IV. Pulse sequence for the implementation of a single 
elementary dissipative map. The necessary operations for the 
reset step are not shown. 



where the two-body terms correspond to interactions be- 
tween spin excitations (or hardcore bosons) located on 
neighboring sites i and i + l. The Hamiltonian dynami- 
cal maps J7i i+i = exp(—i(f>Hi) acting on spins i and i+l 
are up to local rotations equivalent to two-spin MS in- 
teractions, exp(— i<j) o~f crf +1 /4). The implementation of 
the competing Hamiltonian dynamical maps can be re- 
alized with two distinct approaches: (i) The elementary 
two-spin Hamiltonian maps can be implemented sequen- 
tially, in analogy to the sequential implementation of the 
elementary dissipative maps, (ii) Alternatively, since the 
elementary Hamiltonian maps mutually commute, they 
can be implemented by a single global unitary operation, 
acting directly on the entire register of system spins. The 
unitary according to a sum of pairwise interactions be- 
tween neighboring spins can for instance be built up from 
MS gates, which involve pairwise interactions between all 
pairs of spins, by means of refocusing techniques |A4| . Al- 
though one could again try to find numerically optimized 
sequences of gates for the implementation of the compos- 
ite Hamiltonian map, we note that such optimization for 
the unitary acting on the entire register of system spins 
would have to be done for each system size, and becomes 
inefficient for increasing system sizes. It is thus advisable 
to exploit the symmetries of the MS interactions, and to 
systematically construct sequence decompositions. 

Specific experimental implementation - For instance, 
in the case of 3 system spins (with open boundary con- 
ditions), this operation can be decomposed into the fol- 
lowing unitaries according to single- and two-body inter- 
actions 



comp 



= exp 



icj) 



I „Z Z I „z z\ 

(a 1 cr 2 + a 2 a 3 ) 



(Dll) 



i$ t z 
< xi. ! —j{tri 



of + of) exp 



i<j) 



As the ancilla qubit is not required to realize the unitary 
map, and must not be entangled with the system spins 
during the operation, it is spectroscopically decoupled 
during the application of the gate sequence. The experi- 
mental sequences used for the implementation of compet- 
ing Hamiltonian maps for 3 system spins (3+1 ions) and 
4 system spins (4+1 ions) are shown in Tables [vj and VI 
respectively. 



3. Engineering Hamiltonian maps for competing 
coherent interactions 

Competing Hamiltonian dynamical maps are realized 
according to the dimensionless many-body spin Hamilto- 
nian (cf. Section | of the main text) 

JV-l N-l 

H = E^ : = E( 1+o ?x i+o <+i)/ 4 ' ( D1 °) 

i=l i=l 



Appendix E: Experimental characterization of an 
elementary dissipative map and noise model 

1. Modeling an imperfect elementary dissipative 
map 

In the following, we will introduce a theoretical model 
of the elementary pumping step in the presence of experi- 
mental noise. A single elementary dissipative map acting 
on two system ions can be ideally described by the two- 
qubit process matrix Xid which can be straightforwardly 
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Number 


Pulse 


1 


5^1.5,1) 


2 


#(0.5,0.0) 


3 


MS(1- 0.25 k, 0.5) 


4 


541.0,1) 


5 


M5(0.25fe,0.5) 


6 


5,(1.0,1) 


7 


#(0.5,1.0) 



of two depolarizing channels each acting on one of the 
system qubits p M> II(p) — c " 1 7 c" - 
matrix xn as 



ZlZ^iZdepip)) with process 



TABLE V. Pulse sequence for the implementation of a com- 
posite Hamiltonian map in a 3-spin system, according to 
Eq. ( |D11[ >. The Hamiltonian strength is controlled by the 
parameter k, with experimentally implemented values k £ 
{1,0.5}, corresponding to € {tt/2,-k/4}. 



Number 


Pulse 




Act on qubit 2 and 3 


1 


i?(0.5,0.5) 


2 


MS(0.5, 1.0) 


3 


#(0.5,1.5) 




Act on qubit 1,2,3,4 


4 


#(0.5,0.5) 


5 


M5(0.25,1.0) 


6 


5,(1.0,3) 


7 


5,(1.0,2) 


8 


M5(0.25,1.0) 


9 


5,(1.5,3) 


10 


5,(1.5,2) 


11 


#(0.5,1.5) 



TABLE VI. Pulse sequence for the implementation of a com- 
posite Hamiltonian map in a 4-spin system, for strong com- 
peting interactions corresponding to an angle <j> — tt/2. 



calculated from the Kraus map ( |D1[ ) with generating op- 
erators as defined in Eq. ( D2 ) |17|. In order to describe 
the implementation of the elementary dissipative map on 
the three system qubits (two "active" and one spectro- 
scopically decoupled ion), we will assume a process con- 
sisting of the modeled elementary map on the two active 
system qubits and an identity process on the third qubit. 
The noise affecting the elementary map is modeled by 
two independent depolarizing channels acting on each of 
the two "active" system qubits where the fully depolar- 
izing channel on a single qubit i can be written as the 
Kraus map [T7] 



The physical noise acting on the register during the indi- 
vidual gates is certainly more complex than this, but the 
effect of noise on the outcome of any complex algorithm 
can be approximated by depolarizing noise regardless of 
the specific characteristics of the noise [T7] . We model a 
noisy elementary dissipative map with the concatenation 



Xdiss(e) = (1- e)Xid + £Xn- 

In the limit e — > 0, the imperfect process Xdiss reduces 
to the ideal two-qubit process Xid, where in the extreme 
opposite limit of e — > 1 depolarizing noise completely 
dominates and overwrites any effect of the desired en- 
gineered dissipative process Xid- One can now adjust 
the parameter e of this model to the obtained data by 
a numerical optimization. For this, we maximized the 
overlap between the expected output state after a single 
elementary map p e with noise strength e with the actual 
measured state p e xp, 

arg £ maxJ r (p cxp , / 9 e ). 

where the Uhlmann fidelity iF[p\, P2) for the comparison 
of two density matrices pi, pi is used |A8] , We find an 
optimum for a noise parameter of e = 0.27. 



2. Implementation and analysis of an elementary 
dissipative map 

Here, we will provide a more detailed analysis of the 
specific implementation of a single dissipative step. Dur- 
ing the realization of an elementary dissipative Kraus 
map gates act on the two ions encoding the system spins 
i and i+l and the ancilla ion, while all other ions are spec- 
troscopically decoupled. To quantitatively characterize 
the implementation of the Kraus map, as realized by the 
gate sequence shown in Table IV (see also Sec. D2), we 
performed a quantum process tomography on the two 
system qubits. A benchmark for the performance is given 
by the process fidelity with the ideal two-qubit process. 
Since the ideal process is not a unitary process, the Choi- 
Jamiolkowsky process fidelity is a suitable measure |A8| . 
We find a mean state fidelity of J- — 68(1)%. 

In order to identify the leading source of imperfec- 
tions in the implementation, a numerical analysis of the 
actual physical system is performed. We performed a 
Monte Carlo simulation of a single elementary map act- 
ing on three ions, including noise originating from: laser 
frequency and intensity fluctuations, magnetic field fluc- 
tuations, imperfect state preparation, motional heating, 
spontaneous decay and crosstalk of the addressed oper- 
ations. The noise parameters are independently mea- 
sured and we find an overlap between the numerically 
predicted and the measured state of T = 97%. In order 
to determine the dominant error source, we performed 
the simulation multiple times where only a single noise 
source affects the evolution. The results are shown in Ta- 
bic IVIII From this we find that the main error source is 
dephasing due to fluctuations in the laser frequency and 
the magnetic field. 
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Error source 


Overlap with ty + 


All 


77 % 


Addressing 


95% 


Dephasing 


84 % 


Intensity fluctuations 


99% 


Spectator modes 


94% 



TABLE VII. Results for the numerical simulation of a quan- 
tum simulation algorithm, in order to identify the dominant 
experimental error source. The simulation is performed mul- 
tiple times with only a single active error source. From the 
results one can infer that dephasing is the dominant error 
source. 



3. Numerical simulation and long-time dynamics 
under dissipative maps 

Based on the study of experimental errors of an ele- 
mentary dissipative Kraus map above, we have realized 
a numerical study of the dissipative dynamics driving a 
mesoscopic spin systems (N = 10) with initially three 
spin excitations present (m = 3) towards the Dickc 
state |£)(3, 10)). The results are shown and discussed 
in Figure 5a of the main text. To take into account 
the imperfections in the implementation of the elemen- 
tary maps, as discussed above, we assumed a depolar- 
izing noise strength of £diss = 0.02 for each of the two 
spins i and i + 1 involved in an elementary dissipative 
map D i i+ i. Since the implementation of an elemen- 
tary Hamiltonian dynamical map Ui y i+\ on two spins is 
less complex than the realization of a dissipative map, 
we assumed a depolarizing noise strength for f/^j+i of 
( coh = e diss/5- The error reduction protocol to stabilize 
the system within the desired too = 3 excitation sub- 
space can in practice also be only implemented with a 
finite accuracy. Here we assumed both for the global op- 
eration required for one spin extraction procedure and 
for one injection procedure that the whole spin register 
is exposed to a global depolarizing noise channel with a 
probability e = 0.02. 

Long-time dynamics and off-diagonal long-range order 
under imperfect dissipative dynamical maps - As dis- 
cussed above and in the main text, under experimen- 
tal imperfections during the pumping into Dicke states, 
the system initially residing in the too excitation sub- 
space suffers population leakage into the subspaces with 
to 7^ mo. In the long time limit, the population of 
the spin system is homogeneously distributed over these 
rn-subspaces, each subspace populated according to its 
micro-canonical weight, i.e. the number of computa- 
tional basis states spanning the corresponding subspace. 
Within each of the subspaces, not only for to = to , 
the dissipative maps are active and pump the population 
within each subspace towards the corresponding Dickc 
state \D(m,N)). As a result, for weak noise and in the 
long time limit, the system is driven into an incoherent 
mixture of Dicke states of different excitation numbers 



m = 0,...,N. 

More quantitatively, a normalized Dicke state with m 
excitations for N spins is given by 



\D(m,N)) = (to!) 



n- 1 



N 



-1/2 



(S+) m \0)® N (El) 



with S + = X^fcLi a k ■ The number of micro-states of to 

( n\ 

excitations on N sites is . Thus the incoherent 

mixture of Dicke states is given by 



1 N 

w a2 



N 



I' = W M m J \D(m,N)){D(m,N)\ (E2) 

m=0 \ / 



^£(4 (*+)-» 



Interestingly, this incoherent mixture of Dicke states is a 
state with off-diagonal order icxfaj) ^ 0, 



T) ((S+S-) - to) , (E3) 



as we will show in the following. In Eq. (E3) we have 
made use of the symmetry of the Dicke states under 
permutations of the spin indices. In order to evalu- 
ate the global expectation value (S + S~), we decompose 
it as (S+S-) = (S~S + ) + (S z ). Both contributions 
are easily obtained as (S z ) = = ^ m ~ N, and 

(S~ S + ) — (m+l)(N — to). The latter expectation value 
can be obtained using the normalization factors for the 
Dicke states \D(m,N)) and \D(m + 1,N)) with m and 
to + 1 excitations, respectively. Adding both contribu- 
tions, we obtain (S + S~) = m(N + 1 — to). This allows 
us to evaluate the expectation value {o~fo~~) with respect 
to a pure Dicke state: 



(fj O-j )\D(m,N)) 



1 



N(N 
to 
N 



1 



m 
N 



((S+S-) - to) 
1 



1 



j_ 

N 



In the thermodynamic limit, N — > oo, m/N = const., we 
have (cr+af) -> m/N(l-m/N). This expression reflects 
the "particle-hole" symmetry and shows that at complete 
filling (to = iV spin excitations) or in the "vacuum" state 
(to = 0) there is no off-diagonal order because spin exci- 
tations or missing spin excitations cannot be delocalized 
over the spin chain. At half filling to = N/2, where the 
number of micro-states is maximal, the effect of dereal- 
ization and thus the off-diagonal order is maximal. The 
finite-size factor 1/(1 — is to be taken into account 

in small or mesoscopic spin systems, and approaches 1 in 
the thermodynamic limit. 

From this, we can now determine the off-diagonal order 
of the system in the incoherent mixture of Dicke states 
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as given by the density matrix of Eq. ( E3 ) 
identities 

AT / 

V 'm = N2 N ~\ 



Using the 




m 2 = N(N +1)2 



N-2 



m=0 



one obtains 



(f 7 ; 0~<i ) mixture 



As expected, this value does not depend on the initial 
number of excitations anymore, since this information is 
completely lost, once the system has diffused into the in- 
coherent mixture of Dicke states. The off-diagonal order 
in the incoherent mixture of Dicke states assumes a value 
which is independent of the system size N, and remains 
finite in the thermodynamic limit. 

In summary, the experimental imperfections, as de- 
scribed by depolarizing noise and resulting in the non- 
conservation of the excitation-number during the simu- 
lation, lead to a strong decrease of the state overlap be- 
tween the asymptotically reached many-body body state 
of Eq. (|E3} and the ideal "target" dark state \D(m Q ,N)). 
However, from a condensed-matter perspective one can 
state that the imperfections are not too harmful to the 
off-diagonal long-range order, measured by the two-spin 



correlation function (cr^cr - ) for 



> 1 as an order 



•parameter. The off-diagonal order is constantly stabilized 
and re-built by virtue of the repeated application of the 
engineered (though imperfectly implemented) dissipative 
dynamical maps. 



In order to maintain the dissipativcly created off- 
diagonal order in the many-spin system, a crucial prop- 
erty of the excitation number measurement is its QND 
nature: the spin excitation number m has to be deter- 
mined in a way that allows one to only obtain information 
about the total number of excitations in the system, but 
no knowledge about the spatial positions of individual 
excitations along the array. This QND measurement can 
be realized by a global unitary map, which acts on the 
entire register of system spins and an ancillary qubit (see 
Fig. 5b of the main text). Such a unitary is constructed 
in a way that it maps the binary information whether 
the register of system spins is in the correct excitation 
number subspace with to = mo (or not) onto the logical 
state |0) (|1)) of the ancillary qubit. It involves the pro- 
jector PrhP onto the subspace of too spin excitations in 
an array of N spins, and reads 



U(N) 
m 



«p(-i|pW(8a?) (Fl) 

p^^H^+a-pW)®^. 



This equation can be understood as follows: The state of 
the ancilla qubit, initially prepared in |1), is flipped by 
the o~q operation if the system spins are in a state with 
exactly to = too excitations, whereas the ancilla qubit is 
left unchanged otherwise. The QND measurement is then 
completed by a measurement of the ancilla qubit in the 
computational basis, providing the desired information 
on whether m = mo or to too, depending on whether 
the ancilla qubit is measured in |0) or |1). 



Construction of the projectors onto excitation 
number subspaces 



Appendix F: Quantum error detection method: 
Post-selective QND scheme 

1. General idea 

The post-selective error detection method is based on 
a quantum non-demolition (QND) measurement of the 
total number to of spin excitations present in the system 
at the end of the sequence of dynamical maps. Sim- 
ulation outcomes, where due to experimental imperfec- 
tions the ideally conserved initial spin excitation number 
too has changed to a final value to ^ too are discarded. 
This leads to an improvement of the simulation accuracy 
for longer sequences of dynamical maps, at the expense 
of an increased number of experimental runs to obtain 
the same measurement statistics. We remark that in a 
large system, as typical for post-selective techniques, this 
method becomes inefficient as the probability of remain- 
ing within the subspace of initial excitation number Too 
becomes exponentially small and thus only a vanishingly 
small number of "successful" runs enter the measurement 
statistics. 



The projector Piv required for the unitary (|F1| can 



be constructed systematically and efficiently for any too 
and N as follows: One starts from the ansatz 



mo 



N 

E 

k=0 



atkS z 



(F2) 



with S z = 2«=i a t ■ This ansatz assures that the pro- 
jector Pma (i) is diagonal in the computational basis, 
(ii) symmetric under any permutation of spin indices, 
and (iii) does not involve higher powers with k > N as 
such terms are already contained in previous terms with 
k < N due to the property (erf) 2 = 1. Since the compu- 
tational basis states are eigenstates of 5*2, 



S z \0...0) 
5,|0...0,1) 



N\0...0), 
(iV-2)|0...0,l), 



syi...i) = -iv|i...i), 
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and the projector per definition fulfills 
PW|0...0) = 0, 



PW|0...01. 



o... 01 ... i .. 



pW|l l) 

± mo I / 



= o, 



its form is uniquely determined by the following coupled 
system of N + 1 linear equations, 



(in N 2 
1 (N-2) (N-2) 2 



N N \( 

(N-2) N 



1 



-N 



(~N) 2 



\ 



\a N J 



1 




with the only non-zero entry in the (N — m a + l)-th row. 
This matrix equation is readily solved, yielding for the 
experimentally relevant case of three system spins with 
one or two spin excitations the projectors 



p(3) _ 
r \ — 



|011)(011| + |101)(101| + |H0)(110| (F3) 
1 



16 
1 



(9 - 9S Z -Si + S 6 Z ) 



(3-K + af +*$) 

I ^r Z ^ Z I ^ Z ^r Z I „ Z „ Z \ 



and 



P. 



(3) _ 



|001)(001| + |010)(010| 

1 

16 
1 



|100)(100| (F4) 



(9 + 9S Z -Si- S A Z ) 



o (3 + K+a 2 z + a|) 

I ^r Z ^ Z I ^ Z ^r Z I „ Z „ Z \ 0„ Z - Z „ Z \ 

-((r^ +(X 1 a 3 +a 2 a a ) - 6a 1 a 2 a dl ) , 

where S z = °f • Note that the projectors are closely 

related and can be transformed into each other by the 
symmetry operation S z — > — S z , which interchanges the 
role of up- and down-spins, or occupied and empty sites 
in the hardcore boson model, respectively. 



3. Experimental implementation of the QND 
measurement 

As by the total spin excitation number m a collective 
property of the entire spin system is measured, the uni- 
tary of Eq. (Fl ) truly is a many-qubit operation: Equa- 
tions (F3) and (F4) show that the projectors contains 1, 



mapping of Eq. ( Fl I involves interaction terms with up to 
4-body Pauli operators. A decomposition of the unitary 
for the QND measurement U^ =1 into experimentally 
available gates, as obtained using the numerical opti- 
mization algorithm, is shown in Table [VTlT] We note that 
since the QND measurement involves a global unitary, a 
numerical optimization has to be done separately for any 
register size and any particular spin excitation number, 
and furthermore becomes inefficient for increasing regis- 
ter sizes. However, the unitary can be implemented effi- 
ciently without resorting to numerical optimization: For 
a system of N spins, Uma will generally be the product 
of unitaries corresponding to many-spin interactions, at 
most (TV + l)-body Pauli operators. These unitaries can 
be decomposed into sequences of available gates following 
the recipes outlined in JS3|. Although the implementa- 
tion of Um„ becomes experimentally demanding for in- 
creasing system sizes N and in general requires more op- 
erations than numerically optimized circuits, the number 
of required operations for the QND mapping operation 
still grows polynomially with the number of system spins. 



Number 


Pulse 


Number 


Pulse 


1 


#(0.5,-0.5) 


11 


#(0.146,-0.895) 


2 


i?(0.5,0.0) 


12 


M 5(0.375, -1.054) 


3 


540.5,3) 


13 


5 Z (0.364, 3) 


4 


Af 5(0.125, 0.0) 


14 


M5(0.75, -1.054) 


5 


#(0.098,1.0) 


15 


#(1.0,0.0) 


6 


S* (1.636, 3) 


16 


541.818,3) 


7 


M5(0.25,0.0) 


17 


#(0.277,-1.054) 


8 


#(0.136,0.5) 


18 


540.152,3) 


9 


5 Z (0.75,3) 


19 


#(0.5,0.895) 


10 


#(0.113,-1.054) 







TABLE VIII. Pulse sequence for the QND post-selective mea- 
surement of the spin excitation number in a system of 3+1 
ions. 

Due to the considerable complexity of the mapping op- 
eration for the QND measurement, with the optimized 



2 and 3-body spin interaction terms, such that the QND 



circuit involving in total 19 operations (see Table VIII I , 
this method for error detection itself can only be imple- 
mented with a certain accuracy and requires a constant 
resource overhead, independently of the number of dy- 
namical maps applied in the simulation. It can be seen 
from the data shown in Figs. 3a and b, and Fig. 4a of 
the main text, that for short sequences such as e.g. only 
a single elementary dissipative map, where the system 
without error detection remains with high probability 
in the desired excitation number subspace, experimental 
imperfections in the QND measurement itself actually in- 
troduce even more errors on the state of the system (red 
data points) than in the case where it is not applied (blue 
data points) . However, for longer sequences of dynamical 
maps, where the population loss out of the initial exci- 
tation number subspace becomes more and more signifi- 
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cant, the application of the QND post-selective method 
becomes effective and enables a more accurate simulation 
of the system dynamics for longer times. 



qubit. However, in contrast to the QND post-selective 
method, the ancillary qubit is not measured after this 
mapping. The unitary for the mapping reads 



Appendix G: Quantum error reduction scheme 
Excitation number stabilization 



1. General idea 

The second error reduction procedure goes beyond the 
error detection method described above: it does not only 
allow one to detect errors, which have changed the ideally 
conserved spin excitation number to during the quantum 
simulation, but performs an active stabilization of the 
register of system spins in the wanted subspace of spin 
excitation number too- In the previous post-selective case 
the ancilla qubit carried the binary information whether 
the system is in the correct subspace or not. Here, in 
contrast, it is necessary to distinguish between at least 
three cases: (i) The excitation number to is correct, and 
thus no error correction process is required; the excita- 
tion number is (ii) too small or (iii) too large, as illus- 
trated in figure 5b, c in the main text. This information 
cannot be stored in a single ancilla qubit with two log- 
ical states anymore. Although it is in principle possi- 
ble to use a higher-dimensional ancillary system, such as 
multiple ancillary qubits, to store the required informa- 
tion, such an approach would require significantly more 
complex detection and correction algorithms. Here, we 
choose an alternative approach, which still allows one to 
perform the excitation number stabilization with a single 
ancilla qubit if the stabilization process is implemented 
as a two-step process: One stabilization step injects a 
single excitation into the register if there are too few ex- 
citations present (to < mo) and a second stabilization 
step removes an excitation from the system if too many 
are present (m > mo), as schematically shown in Fig. 5c 
of the main text. Similarly to repetitive quantum error 
correction, the state of the ancilla qubit has to be reset 
in between the two steps. As the underlying construc- 
tion of the excitation injection and the extraction step 
is very similar, we focus in the following on a detailed 
description of the protocol for the excitation extraction 
procedure. 

Conceptually, the injection protocol consists of two 
main parts: First, similar to the QND measurement 
above, the binary information whether or not there are 
too many spin excitations in the system is coherently 
mapped onto the two logical states of an ancillary qubit. 
Second, conditional on the state of the ancilla qubit, 
which is acting as a quantum controller, a feedback pro- 
cedure is applied to the system spins. This quantum 
feedback procedure extracts one of the (possibly multi- 
ple) superfluous spin excitations and stops once this has 
been achieved. 

The first part is realized by a unitary mapping act- 
ing on the entire register of system spins and the ancilla 



U, 



(JV) _ 



JV 



m>mo 



exp —i- 



E 3 

v j=m + l 



(N) 



(Gl) 



and can be understood as follows: The state of the an- 
cilla qubit prepared in |1) is flipped if and only if there 
are too many spin excitations present in the system; here 
the system operator J2f= mo +i denotes the projector 
onto the subspace of states containing strictly more exci- 
tations than the desired value Too- It is the sum over all 
projectors Pj N ^ onto subspaces with j > mo excitations 
in a system of N spins, and each of these projectors can 
be readily constructed following the procedure described 
in Sec. ED 

After the mapping operation between the system spins 
and the ancilla, the second part of the protocol deals with 
the actual extraction of a spin excitation from the sys- 
tem. This step faces the difficulty that due to the QND 
character of the first mapping step, the state of the an- 
cilla qubit only encodes information about whether too 
many excitations are present in the system, but not on 
which sites the excitations are located. In order to mini- 
mally disturb the state of the system spins it is desirable 
to devise an extraction scheme with the following proper- 
ties: (i) An excitation is extracted in a minimally invasive 
way, i.e. under an extraction of an excitation from a cer- 
tain site, off-diagonal order among system spins of the 
rest of the chain is maintained as far as possible, (ii) 
Second, the scheme should effectively hold and not fur- 
ther alter the state of the system spins once an excitation 
has been successfully extracted. We have developed and 
implemented a scheme, which satisfies both criteria, and 
which exploits a combination of spectroscopic decoupling 
and optical pumping of the ancillary qubit. 



2. Steps of the excitation extraction procedure 

Let us now describe in more detail the steps of the 
excitation extraction procedure, which is illustrated in 



Fig. A4 for a system of three spins. In step (1) the map- 
ping introduced above is applied, so that the information 
whether a correctable error (i.e. too many spin excita- 
tions in the system) is present (or not), is encoded in the 
ancilla state |0) (|1)). If no correctable error is detected 
(to < mo) then the subsequent operations should not 
affect the system qubits. This can be realized by remov- 
ing the ancilla qubit incoherently from the computational 
subspace (step (2)) as will be discussed below. In step (3) 
an extraction attempt on the first site of the spin regis- 
ter is performed by swapping the ancilla, which is now in 
the state |0) corresponding to an empty site (no spin ex- 
citation), with the first qubit of the system. In terms of 
spins this operation corresponds to a swap ("flip-flop") 



21 



Excitation number stabilization procedure 
(spin excitation / hardcore boson extraction protocol) 



ancilla 



system 



(1) QND map: m > 1 ? (1) 




1 2 3 

ideally: m = 1 excitation 



if m = 2 or m = 3 



2 3 m = 0: error 

(not correctable now) 



m = 2: error 
I — 9_i L*J (1 exc. too much) 

o or 1 2 3 

•JL«JL*J m = 3: error 
1 2 3 (2 exes, too much) 



(2) removal 
pumping 
state |2> | 



of the ancilla by 
to'parking" 




no further 
effect of 
steps 
(3) -(5) 
on the 
system 



^ final state: 

LlJL 



no error 

or uncorrectable 
m=0 error 



12) 

(3) exc. extraction 

intent from site #1 / 1 exc. 

extracted 




(4) removal of the ancilla 
to "parking" 
state |2) 





L»f-tJ 



no further 
effect of 
step (5) 



(5) exc. extraction 
from site #2 



(5) 



1 exc. extracted 
with certainty 



(6) reset of the 
ancilla to |0) 



Suj final state: 



or 



1 



2 3 

II I 

2 3 



J I error completely corrected 

exc. number reduced from 
2 3 m=3 to m=2 




1 2 3 

error completely 
corrected 



Figure A4. Logical tree illustrating the complete excitation extraction procedure for a 3-spin system with ideally a single 
excitation present, and for various (unknown) initial states of the system spins. Blue circles correspond to spin excitations 
(hardcore bosons) localized on particular sites, whereas ellipses indicate situations where excitations (hardcore bosons) can be 
delocalized over several lattice sites. The error reduction takes place in six steps: By a QND map (1) the information whether 
or not too many excitations are present in the system, is mapped onto an ancilla qubit. (2) In case of the ideal excitation 
number or too few excitations, the ancilla is removed from the computational space and the protocol effectively halts. In the 
opposite case, a swap operation (3) is applied to try to extract an excitation from system site # 1. If successful, after the 
subsequent removal of the ancilla (4) the protocol halts. Otherwise, (5) an extraction by a swap of the ancilla with system site 
# 2 extracts an excitation. (6) Finally, the ancilla is reset to its initial state |0) for subsequent error reduction steps. See also 
text for a more detailed explanation. 
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process, under which |0)o|l)i — > |l)o|0)i, i.e. the spin 
excitation is coherently exchanged between the ancilla 
qubit and the first system spin. In the language of hard- 
core bosons, the process can be interpreted as a tunneling 
of a boson from the site #1 onto the previously empty 
ancilla site (see step (3) in Fig. A4). If after this process 
the ancilla qubit is in the state corresponding to an oc- 
cupied site (spin excitation present), a single excitation 
has been successfully extracted from the system and the 
process should halt. As will be shown below, this can 
be performed by removing the population of the ancilla 
qubit from the computational subspace (step (4)). How- 
ever, if the ancilla is still in the state |0) corresponding 
to an empty site after the swap operation, then the first 
system qubit was in the state corresponding to an un- 
occupied site and thus no spin excitation (or hardcore 
boson) could extracted, |0)o|0)i — ► 1 0) q 1 0) i . Thus, this 
procedure is repeated (step (5)) until the ancilla qubit is 
found in the state |1) corresponding to an occupied site. 
In the case of three sites at most two extraction rounds 
are required (steps (3), (4), and (5)), before the ancilla 
is reset to |0) (step (6)) for subsequent rounds of error 
reduction. 



Open vs. closed-loop error reduction scheme — The de- 
scribed extraction (and injection) protocols are realized 
in an open-loop fashion, i.e. the ancilla remains unob- 
served during the whole procedure. This comes at the 
cost that all pulses (for steps (1) - (6) in the above exam- 
ple) have to be physically applied in every run of the error 
reduction procedure, even if the ancilla has been already 
be removed from the computational subspace and in prin- 
ciple no further operations would be required. However, 
the described protocol could easily be modified into a 
closed-loop control scenario: In this case one would per- 
form individual measurements of the ancilla qubit after 
the QND map (1) as well as after the coherent swap oper- 
ations (steps (3) and (5)). The outcomes of these ancilla 
measurements would then yield the classical information 
on whether no error correction is required after step (1) 
or an excitation has already been extracted successfully 
in steps (3) or (5). This information could then be used 
for classical feedback on the quantum system in the sense 
that one can externally decide whether further steps are 
required in the error reduction protocol, or whether one 
can stop the current run of the error reduction protocol 
and no further operations have to be applied. An advan- 
tage of this modified scheme is that dissipative removal 
operations of the ancilla qubit out of the computational 
subspace are not required. Experimentally, the price to 
pay in this modified scheme would be a number of fluo- 
rescence measurements, which are slow compared to the 
application of coherent gates, and the requirement to re- 
cool the relevant external degrees of freedom of the ion 
chain after each measurement as recently demonstrated 
in [AT] , 



3. Experimental implementation 

In the experiment, the error reduction protocol for a 
system of three spins, using 3+1 ions was implemented. 
We considered the case of having ideally m = 1 ex- 
citation present on the three sites. For the extraction 



(3) 

step, the four-qubit unitary (cf. Eq. (Gl|) to de- 

tect whether too many excitations (i.e. m = 2,3) are 
present in the system or not, has been obtained by nu- 
merical optimization. The resulting circuit consisting of 
15 gates is listed in Table |IX| On the other hand, for the 
injection part, the unitary 



U, 



(AT) 



m<mo 



exp 




(G2) 



which interrogates the system whether too few excita- 
tions are present, is required. For the case mg = 1 and 

(3) 

N = 3, the unitary U^i realizes a spin flip of the ancilla 
qubit only if the three system spins are in the state 1 000) . 
This operation is equivalent to a triple controlled-NOT 
operation with the three system spins playing the role of 
the control qubits and the ancilla the target qubit. As 
we could not obtain a satisfying circuit decomposition 
for this unitary by means of the usual numerical opti- 
mization algorithm, we did not try do directly optimize 

(3) 

^m<i> but allowed the unitary to add arbitrary phases 
to states lying outside the desired excitation number sub- 
space m = 1. Under these relaxed conditions, the nu- 
merical optimizer delivered a circuit decomposition of 19 
operations as shown in Table |X| 



Number 


Pulse 


Number 


Pulse 


1 


#(0.25,0.5) 


9 


5,(1.0,3) 


2 


5,(1.0,3) 


10 


#(0.25,0.5) 


3 


#(0.25,0.5) 


11 


MS(0.5, 0.5) 


4 


M5(0.25,0.5) 


12 


#(1.5,0.5) 


5 


#(1.75,0.0) 


13 


#(0.125,0.5) 


6 


5,(1.0,3) 


14 


5,(1.0,3) 


7 


#(0.25,0.0) 


15 


#(1.875,0.5) 


8 


#(1.75,0.5) 







TABLE IX. P ulse sequence to implement the QND unitary 



^m>i ( CL Eq. (Gil) as part of the excitation extraction pro- 



tocol. It maps the information whether there are more than 
m — 1 spin excitations present in a system of three spins, 
onto the ancilla qubit. 

Implementation of the actual excitation extraction - 
We will now discuss the physical mechanism underlying 
the implementation of the removal of a spin excitation 
(or hardcore boson), as described above and illustrated 
in Fig. |A5| This excitation removal step takes place af- 
ter the QND mapping (1) and a decoupling operation 
on the ancilla qubit (2) described below. As mentioned 
above, the extraction step (3) relies on a swap operation 
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Number 


Pulse 


Number 


Pulse 


1 


#(1.75,0.5) 


12 


#(0.375,0.0) 


2 


&(1.0,3) 


13 


^(1.0,3) 


3 


#(1.75,0.5) 


14 


#(1.625,0.0) 


4 


M5(0.25,0.0) 


15 


S*(0.5,3) 


5 


#(1.875,0.0) 


16 


M 5(0.25, 0.0) 


6 


5 Z (1.0,3) 


17 


#(0.125,0.0) 


7 


#(0.125,0.0) 


18 


S*(1.0,3) 


8 


#(0.25,0.5) 


19 


#(1.875,0.0) 


9 


5 Z (1.0,3) 


20 


#(0.125,0.5) 


10 


#(1.75,0.5) 


21 


S*(1.0,3) 


11 


MS(0.25, 0.0) 


22 


#(1.875,0.5) 



TABLE X. P ulse sequence to implement the QND unitary 



Um<i (°f- Eq. (G2 i) as part of the excitation injection proto- 



col. It maps the information whether there are m = instead 
of ideally m = 1 spin excitations present in a system of three 
spins, onto the ancilla qubit. 




• - 



- halt? J — > 


halt? 


















■ 








Figure A5. a, Schematic view of the error detection and 
excitation removal process (cf. Fig. 5c of the main text), b, 
The excitation removal process on 3 system qubits can be 
performed by swap operations and effective halting condi- 
tions realized by a dissipative decoupling of the ancilla qubit. 
c, Swap operations (G3l between the ancilla qubit and one 



system spin are implemented by two effective two-qubit MS 
operations, which are build up from 4 MS operations acting 
on three ions, interspersed with refocusing operations. 



Dissipative decoupling of the ancilla ion 
— 4 2 P. 




^4 2 P 




4 2 P, 



3 2 D 

|o> 



5/2 



12) _ 




-1/2 
-3/2 
-5/2 



4 2 S, 



1/2 
-1/2 



Figure A6. The dissipative decoupling process of the ancilla 
qubit used in the injection/removal processes. The goal is 
to uni-directionally remove population in the computational 
state |0) (red triangular marker) and add it to electronic pop- 
ulation that is possibly already present in the "parking" state 
1 2} (green rectangular marker). Population in the computa- 
tional state |1) (blue circle) should be left unchanged. The 
sequence: (i) First, the population in states |1) and |2) is co- 
herently swapped by a 7r-operation. (ii) Population in |0) is 
coherently transferred to the state 4Si/2(rnj = +1/2). (iii) 
Optical pumping from this state towards the |1) state. In 
this step, the populations that were at the beginning of the 
sequence present in \2) and |0) are added up and temporarily 
accumulate in state |1). (iv) Finally the populations in states 
|1) and 1 2} are again swapped coherently. 



exchanging the excitation of the ancilla qubit with one 
of the system spins, say #1. This swap operation is re- 
alized by applying a 7r-pulse, according to the "flip-flop" 



Hamiltonian 



H, 



01 



:( cr o cr i 
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Number 


Pulse 


1 


M5(0.25,0) 


2 


s.(i.o,i) 


3 


MS'(0.25,0) 


4 


MS(0.25, 0.5) 


5 


&(i.o,i) 


6 


MS(0.25,0.5) 



TABLE XI. Pulse sequence for a swap operation used for ex- 
citation injection or removal, see Fig. [A5ji. Ion j is the spec- 
tator ion, which is not spectroscopically decoupled, but is not 
intended to participate in the swap operation. 



The resulting unitary 



U i =exp 



(G3) 



corresponds to the application of tw o fully-entangling x- 
and y-type MS gates (see Eq. ( |C3[ )) applied to the an- 
cilla ion and the first system ion. This two-ion MS gate 
can be realized by the global bi-chromatic light field, if 
all ions but the ancilla qubit and the system qubit #1 
are spectroscopically decoupled. However, since step (3) 
involves the ancilla qubit and system ion #1, and step 
(5) the ancilla qubit and system ion #2, it is from an 
experimental point of view more convenient to spectro- 
scopically decouple only ion #3, but keep both ions #1 
and #2 in the "active" qubit states during the extraction 
procedure. Doing so, each of th e x— and y-type MS-gates 
required for the operation ( G3 ) acting on the ancilla (in- 



dex #0) and system ion #1 can be realized by two MS 
gates acting on the three ions, interspersed by a 7r-pulse 
5z(l,2) applied to the second system ion, as shown in 
Fig. |A5| The application of the "refocusing" pulse in 
between the MS gates leads to an effective decoupling of 
ion #2 that is not supposed to participate in the unitary 
(G3| |A4| . This technique is also applied to realize the 
swap operation required in step (5), which acts on the 
ancilla qubit and the second system. Using this decou- 
pling technique based on refocusing pulses, no additional 
spectroscopic decoupling operations are required between 
the two swap operations corresponding to steps (3) and 
(5). The experimentally employed sequence of 6 gates is 
listed in Table ED 

In Sec. IG 21 we have outlined that the ancilla will be 
removed from the computational subspace in step (2) if 
no error is to be corrected in the current error reduc- 
tion round. In this case, if the ancilla has been removed 
from the qubit subspace, the two-qubit MS interactions 
appearing in (G3) act only on a single qubit. However, 
in this "pathological" case where the bi-chromatic laser 
fields used for the generation of the effective spin-spin 
interactions of the entangling MS gate of Eq. ( C3 1 , are 



[2] , As a consequence, in this case - despite the 
fact that the MS gate laser pulses of steps (3) and (5) 
are physically applied to the ions - they do not alter the 
state of the systems spins, as desired. 

Dissipative decoupling of the ancillary qubit - As ex- 
plained above, the ancilla qubit shall be removed from 
the computational subspace if either (i) no correctable 
error is detected in the current round of error reduction, 
or (ii) a superfluous spin excitation has been successfully 
extracted from the spin system. In these cases, the re- 
moval of the ancilla from the qubit subspace guarantees 
that the error reduction protocol effectively halts (see the 
logical tree shown in Fig. A4 1 and the state of the sys- 



tem spins is no further modified, even though after the 
removal of the ancilla MS gate laser pulses are applied to 
the ion string. In steps (2) and (4) of the spin extraction 
protocol, the ancilla qubit is removed from the compu- 
tational state |1) into the "parking" state |2), which we 
encode in the ZD^^ijnj — —5/2) electronic state (see 
Fig. A6). This removal into |2) takes place if the ancilla 
ion resides at these instances in |1), if it is in |0) its state 
remains unaffected. 

For this removal it is not possible to use the coher- 
ent spectroscopic decoupling technique as used for the 
implementation of the elementary Kraus maps. The rea- 
son is that this would lead to errors in the protocol: for 
instance, imagine no correctable error is present in the 
system and thus the ancilla is brought to the decoupling 
state 1 2) in step (2). If then another ancilla removal 
were performed coherently (step (4)), the ancilla would 
be transferred back from |2) to the computational state 
|1). This is clearly unwanted as in this case with the an- 
cilla returned to |1), the subsequent swap operation (5) 
would be performed by mistake. This unwanted behavior 
can be avoided if the transfer of the ancilla from |1) to 
1 2) is realized dissipatively, by an optical pumping pro- 
cess, which bares similarities with the incoherent reset 
of the ancilla qubit for the elementary dissipative maps. 
Such uni-directional, incoherent pumping process from 
|1) to 1 2) guarantees that once the ancilla has reached the 
"parking" state |2), it will in subsequent steps never re- 
turn to the computational subspace. This optical pump- 
ing process for a removal process of the ancilla qubit from 
one of the computational basis states into |2) is illustrated 
and described in more detail in Fig. |A6| 



applied to a single ion only, these realize up to negligi- 
ble corrections the identity operation on this single ion 



In the spin excitation protocol the ancilla is removed 
from |1) (corresponding to an occupied site) to |2), 
whereas in the spin injection protocol it has to removed 
from |0) (empty site) to |2). For both scenarios, the 
pulse sequence outlined in Fig. |A6| can be employed; 
a 7r-pulse resonant with the qubit transition of the an- 
cilla, applied before and after the dissipative decoupling 
sequence, exchanges the roles of the two computational 
states |0) -O- |1) and thereby allows one to switch between 
the spin extraction and injection scenario. 
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4. Experimental results for the stabilization 

The active excitation number stabilization procedure 
can be best tested when applying it to a state that has a 
considerable amount of the population outside the sub- 
space with the correct excitation number uiq. We ap- 
plied Hadamard operations on the three system qubits 
initially in |000) to prepare the initial state \ipo) — 
l/v / 8(|0} + ll))® 3 , which is an equal-weight superposition 
of all eight three-qubit computational basis states, each 
occurring with probability 1/8. The measured and ideal 
density matrix of this initial state is shown in Fig. \A8\ i. 
For m = 1 the states |001), |010), |100) span the sub- 
space with the desired spin excitation number. Thus, 
the initial fraction of population in this subspace is 3/8 
as shown in Fig. 5d of the main text. We then performed 
the excitation extraction protocol according to the proto- 
col outlined above and summarized in Fig. |A4| Ideally 
this protocol extracts one excitation from the compo- 
nents of the initial state, which contain two or three spin 
excitations, m = 2 — > m = 1 and m = 3 m = 2, 
thereby pumping the population corresponding to these 
states into the subspaces with one excitation less. It is a 
crucial property of the error reduction protocol that the 
coherences within the subspace of the "correct" excita- 
tion number m — 1 are ideally preserved, as the dynamics 
within the desired simulation subspace should be affected 
as little as possible. The component 1 000) of the initial 
state corresponds to a state with zero, i.e. too few spin 
excitations, and thus to an error which is not corrected 
by the spin excitation procedure. 

In Fig. 5d in the main text, the ideal and the mea- 
sured populations in all four excitation number subspaces 
at the end of the extraction protocol are shown. From 
this information one can deduce that the protocol within 
experimental accuracy realizes the pumping of popula- 
tions between the different excitation number subspaces 
as expected. To infer whether the coherences are pre- 
served, we measured the three-qubit density matrices af- 
ter the first excitation extraction attempt from the first 
site (step (3) in Fig. A4) and after the second extraction 
attempt from the second site (step (5)). The measured 
and the ideal density matrices are shown in Fig. |A8b 
and c, where the relevant coherences within the uiq = 1 
subspace are highlighted in red color. From the com- 
parison of measured and ideal density matrices it can be 
seen that the coherences in the mo excitation subspace 
are well-preserved. 

The complete error correction protocol consists of the 
excitation extraction and injection procedures. For the 
same initial state \"ipo), we implemented the injection pro- 
cedure, which ideally only acts on the 1 000) state, and 
pumps the population from the m = into the m = 1 
subspace. Again, the coherences within the m = 1 sub- 
space should be preserved under this procedure, yielding 
as a result of the injection protocol the ideal density ma- 
trix shown in the right part of Fig. |A8[ 1. Comparison 
with the measured density matrix shows that most of the 



population is pumped out of the m = subspace, and 
that the initially present coherences within the m = 1 
subspace are reasonably well conserved. 

The ultimate goal of any error reduction protocol is 
certainly to increase the performance of a complex al- 
gorithm. As a step in this direction, we integrated the 
excitation removal protocol into the simulation sequence 
for dynamics according to composite dissipative maps 
with 3+1 ions. In Fig. |A7| we compare the probabilities 
for all excitation numbers m when (i) no error reduction 
technique is used (blue data points), (ii) with the post- 
selective QND measurement applied (red), and (iii) with 
the excitation removal procedure included in the simula- 
tion (green). It can be seen that the removal procedure 
has a higher overhead due to its considerable complex- 
ity. Nevertheless, as an indication of its usefulness, a 
slower decay of the probability of finding the population 
in the desired subspace for m = 1 can be observed, for 
the case in which the stabilization procedure is applied, 
compared to the case without any error correction. This 
indicates that the stabilization procedure indeed works 
qualitatively as intended, when it is incorporated into 
the actual simulation. 
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Figure A7. Experimental study of the error reduction proto- 
col, incorporated into the actual quantum simulation of com- 
posite dynamical maps with 3+1 ions, with an initial and ideal 
excitation number of m = 1. The four plots show the fraction 
of the population in the subspaces with zero, one, two and 
three spin excitations, as a function of the number of elemen- 
tary dissipative dynamical maps applied. For comparison, 
blue data points corresponds to the population under dissi- 
pative maps, without error detection or reduction technique 
applied. Red data points correspond to the case where after 
the final elementary dissipative map the QND post-selective 
measurement was applied. Green data points corresponds to 
the case where the extraction protocol has been applied at the 
end of the simulation. Bars correspond to the theory, where 
depolarizing noise in the elementary dissipative maps is taken 
into account (cf. Appendix Ell. 



2G 



initial state - experiment 



initial state -theory 



Extraction attempt site 1 
experiment 



Extraction attempt site 1 
theory 







Extraction attempt site 1+2 
experiment 



Extraction attempt site 1+2 
theory 



injection attempt site 1 
d experiment 



injection attempt site 1 
theory 





DDD 

DDS 

DSD DSS ^-^T A/ DSD 
SDD V*/,.,, 

SSD DDD 
SSS uuv 





Figure A8. Reconstructed (left) and ideal (right) density matrices of the removal and injection process in a 3+1 ion experiment. 
Populations and coherences within the desired excitation number subspace with m = 1 are high-lighted by red bars, a, The 
initial state is an equal-weight superposition of all eight computational basis states, b, The state after the system after a 
spin excitation removal attempt from the first site, c, The state of the system after the second swap operation to remove an 
excitation from the second site. After both steps, the coherences shown in red are well-preserved, d, The state of the system 
after a spin excitation injection attempt, starting again in the equal-weight superposition state shown in a. The population 
from the \DDD) = |000) state is depleted, and coherences within the m = 1 subspace are reasonably well preserved. 
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Figure A9. Experimental results of dissipatively 
induced derealization through composite dynamical 
maps with 3+1 ions. The results from an ideal model 
are shown in light-blue bars whereas those from a model in- 
cluding depolarization noise are indicated by dark-grey bars. 
Blue rectangles indicate the experimentally observed dynam- 
ics without any correction scheme whereas red rectangles in- 
clude a post-selective error detection scheme (error bars, la). 
Overlap fidelity purity, population in the m — 2 subspace, 
and off-diagonal order in a 4-spin quantum simulation with 
4+1 ions, studying purely dissipative dynamics that induces 
pumping towards the Dicke state \D(2, 4)} as shown in figure 
3 in the main text. 
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Figure All. Experimental results for competing dis- 
sipative and coherent dynamics with 3+1 ions with 
a single excitation. The results from an ideal model are 
shown in light-blue bars whereas those from a model in- 
cluding depolarization noise are indicated by dark-grey bars. 
Blue rectangles indicate the experimentally observed dynam- 
ics without any correction scheme whereas red rectangles in- 
clude a post-selective error detection scheme (error bars, la). 
Data from simulated dynamics with Hamiltonian competi- 
tion for 3 spins (3+1 ions), but only a single spin excitation 
present. As physically expected, the data confirms that in 
this case the Hamiltonian dynamical map has no effect. 
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Figure A10. Experimental results for competing dis- 
sipative and coherent dynamics with 3+1 ions The re- 
sults from an ideal model are shown in light-blue bars whereas 
those from a model including depolarization noise are indi- 
cated by dark-grey bars. Blue rectangles indicate the exper- 
imentally observed dynamics without any correction scheme 
whereas red rectangles include a post-selective error detection 
scheme (error bars, la). Overlap fidelity, purity, population in 
the m — 2 subspace, and off-diagonal order in a 3-spin quan- 
tum simulation with 3+1 ions. The dynamics corresponds 
to dissipative maps and coherent Hamiltonian competition as 
shown in figure 4 in the main text. 



Figure A12. Experimental results for competing dis- 
sipative and coherent dynamics with 3+1 ions with 
weak interaction The results from an ideal model are shown 
in light-blue bars whereas those from a model including depo- 
larization noise are indicated by dark-grey bars. Blue rectan- 
gles indicate the experimentally observed dynamics without 
any correction scheme whereas red rectangles include a post- 
selective error detection scheme (error bars, la). Data from 
simulated dynamics with weaker competing Hamiltonian dy- 
namics (<f) — 7r/4) in a 3-spin system with two spin excitations 
present in the system (3+1 ions). 
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Appendix H: Additional experiments and data 
analysis for competing dissipative and coherent 
dynamics 

In Figs. 3 and 4 of the main text, the experimental 
results for systems of 3+1 and 4+1 ions are depicted. 
The attentive reader noticed that in Fig. 4a in the main 
text the data point for 5 maps is missing. This is due 
to the fact that in view of the length of the algorithm, 
the memory of the current experimental control system is 
not sufficient to generate the required pulse sequence for 
4 elementary dissipative maps and one composite Hamil- 
tonian dynamical map. It is however possible to generate 
the sequence for two elementary maps and one Hamilto- 
nian dynamical map, and to repeat this sequence twice 
(data point for 6 maps). 

Here, we extend the experimental analysis that was 
omitted from the main text by showing different mea- 
sures for the already presented data and also additional 
datasets: In Figs. |A10 and A9 we add the measures 
purity and off-diagonal order for the data already pre- 
sented in Figs. 3 and 4 of the main text. The purity 
Trp 2 is a measure for how close the measured state is 



to a pure state. Off-diagonal order measures the coher- 
ences between neighboring sites as the expectation value 
of the operator V ■ evaluated within the mo- 

excitation subspace. This parameter emphasizes the ef- 
fect of the competing Hamiltonian dynamics as it changes 
the sign from positive to negative after the application of 
a Hamiltonian map, as shown in figure [A10| 

Figure [A 1 1 1 shows an alternative dataset which demon- 
strates that the coherent competing Hamiltonian dynam- 
ical map in a system of three sites has - as physically ex- 
pected - no effect if only a single excitation is present in 
the system. Comparison with the data of Fig. 4a in the 
main text, underlies the significance of the decrease in 
the overlap with the overlap with the target Dicke state 
|D(2, 3)), and this effect is indeed caused by the compet- 
ing Hamiltonian dynamics. 



Figure [A 1 2 1 shows a dataset with competition and two 
spin excitations present in the system. Compared to the 
analysis in the main text, the competition strength is now 
set to 4> = 7r/4 instead of ir/2. As expected, this leads to 
a reduced effect of the competing Hamiltonian maps on 
the dissipatively created order. 
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